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1  .  I  r;f  rnciiii' t  ion 

in-.piovod  capabilities  to  specity  the  moisture  field  (including  water 
vapor,  cloud,  and  precipitation)  are  potentially  afforded  by  an  optimally 
designed  combination  of  satellite-based  sensors.  The  implications  for 
numerical  weather  prediction  (NWP)  are  promising,  since,  radiatively,  water  in 
its  manv  forms  is  the  most  active  constituent  of  the  troposphere  (Isaacs  et 
al.,  1986).  While  infrared  cloud  contamination  may  often  be  circumvented  by 
exploiting  the  high  horizontal  resolution  of  infrared  sensors,  microwave 
temperature  retrievals  suffer  little  degradation  from  the  presence  of  cloud. 
This  advantage  also  applies  to  some  extent  to  the  retrieval  of  the  atmospheric 
moisture  profile  from  millimeter  wave  data. 

A  strong  feature  with  sufficient  line  strength  for  this  application,  is 
the  rotational  line  of  water  vapor,  located  at  183.31  GHz.  Two 

ins t riunonts  liavc  been  proposed  to  utilize  the  183  GHz  absorption  line  to 
obtain,  water  vapor  profiles  in  the  1990' s  time  frame.  These  are;  (a)  the 
civilian  Advanced  Microwave  Sounding  Unit  (AMSU)-B  package,  a  five  channel 
radiometer  to  be  flown  as  part  of  the  first  NOAA-NEXT  satellite,  and  (b)  the 
Defense  Meteorological  Satellite  Program  (DMSP)  SSM/T-2  millimeter  wave 
moisture  sounder,  an  enhancement  of  the  current  operational  SSM/T-1  microwave 
temperature  sounder. 

This  study  presents  the  results  of  state  of  the  art  observing  system 
simulation  experiments  (OSSEs)  designed  to  assess  the  impact  of  of  the  Special 
Sensor  Microwave  183  GHz  (SSM/T-2)  water  vapor  sounder  and  the  associated 
SSM/T-1  temperature  sounder  on  NWP.  A  novel  feature  of  the  present  study  is 
the  comprehensive  approach  we  have  taken  to  simulating  the  SSM/T  microwave 
(T-1)  and  millimeter  wave  (T-2)  data.  Our  study  is  based  on  the  nature  run 
prepared  by  the  ECMWF  and  simulated  FGGE  data  base  prepared  by  NMC  as  de¬ 
scribed  by  Dey  et  al,  (1985)  for  the  period  10-30  November  1979.  The  forecast 
analysis  system  for  the  OSSE  is  the  AFGL  Global  Spectral  Model  (GSM)  (Brenner 
et  al.,  1984)  and  the  AFGL  Statistical  Analysis  Program  (ASAP)  (Norquist, 

1986).  Four  observing  system  configurations  are  investigated:  (1)  in  experi¬ 
ment  NOSAT  only  the  conventional  upper  air  soundings  and  aircraft  reports  are 
used;  (2)  in  experiment  STATSAT,  the  currently  operational  civilian  tempera¬ 
ture  soundings  (TOVS)  are  added;  while  (3)  in  experiment  SSMSAT  the  SSM/T 
temperature  and  moisture  soundings  are  added  to  the  NOSAT  configuration; 
finally  (4)  in  experiment  SSM-tTOVS,  both  the  SSM  and  TOVS  data  are  added. 
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As  an  aid  to  interpretation,  our  OSSE  results  are  calibrated  against 
results  from  (real  data)  observing  system  experiments  (OSEs)  described  in  more 
detail  by  Louis  et  al .  (1988).  In  addition  to  conventional  assessment  indices 
(such  as  geopotent j ai  height  fields,  etc.),  our  analysis  of  the  results  also 
focuses  on  moisture  related  variables.  We  Identified  two  measures  of  impact 
of  particular  interest:  relative  humidity  and  cloud  amount.  Due  to  the 
importance  of  moisture  fields  to  the  prediction  of  electro-optical  systems 
effectiveness,  it  is  useful  to  study  a  variable  directly  related  to  water 
vapor.  Relative  humidity  is  the  parameter  most  commonly  used  as  an  input  to 
propagation  models  (Moore  and  Peterson,  1984).  This  assessment  will  be 
important  for  situations  with  moderate  to  high  relative  humidities. 

For  situations  near  saturation,  cloud  prediction  indices  are  most 
important.  Therefore  we  have  examined  the  rms  and  bias  error  of  cloud 
fraction.  Since  cloud  amounts  were  not  saved  during  the  nature  run,  truth  is 
based  on  the  diagnosis  of  large  scale  (i.e.  non  convective)  cloudiness  from 
the  nature  run  relative  humidity  fields.  We  find  that  the  statistical 
measures  of  impact  based  on  cloud  amount  are  closely  related  to  those  based  on 
relative  humidity.  We  note  Chat  special  attention  is  given  to  the  simulated 
retrieval  errors  in  the  presence  of  cloud.  Cloud  and  retrievals  ove^r  land  are 
two  situations  where  problems  with  Che  183  GHz  retrievals  are  expected  based 
on  previous  retrieval  simulation  studies  (Isaacs  and  Deblonde,  1985). 

The  plan  of  this  report  is  the  following:  Section  2  discusses  OSSE 
methodology  in  general  and  develops  the  methodology  for  realistic  simulation 
of  the  SSM  sensors.  Background  material  describing  the  scientific  basis  for 
the  183  GHz  sensor  have  already  been  reviewed  by  Isaacs  (1987).  This  material 
is  not  repeated  here.  Section  3  describes  the  reference  atmosphere  and 
simulated  conventional  data.  Section  4  describes  the  sensor  simulation  study 
of  the  SSM/T  statistical  (i.e.  D-raatrix)  retrievals,  and  our  plan  for 
simulating  SSM/T  data  for  the  OSSE.  The  simulated  SSM/T  data  observational 
errors  are  based  on  the  results  of  our  sensor  simulation  retrieval  study. 

These  results  are  examined  in  Section  5.  Section  6  contains  the  specific  OSSE 
design  we  employed.  Results  are  reported  in  Section  7.  Finally,  Section  8 
contains  a  summary  and  our  main  conclusions. 

The  experiments  described  here  were  conducted  in  tandem  with  an  OSSE  to 
determine  the  impact  of  Doppler  wind  lidar  Instrument  or  WINDSAT.  The  WINDSAT 


study  makes  use  of  the  same  baseline  OSSEs  and  OSEs  as  the  current  study.  The 
WINDSAT  study  is  described  in  a  companion  report  (Scientific  Report  No.  7)  of 
the  current  effort.  Each  report  may  be  read  independently:  To  achieve  this 
we  have  repeated  some  common  material  in  both  reports. 

2 .  Background 

OSSE  studies  run  the  gamut  from  simple  insertion  of  grid  point  values 
into  unsophisticated  models  to  the  use  of  complex  radiative  transfer  models  to 
simulate  data  for  operational  forecast  analysis  systems.  The  current  study  is 
relatively  sophisticated  and  state  of  the  art.  As  a  reference  point  we 
discuss  in  Section  2.1,  considerations  which  might  lead  to  an  ideal  OSSE.  We 
then  describe  our  methodology  for  simulating  remotely  sensed  data  (Sec¬ 
tion  2.2).  This  methodology  is  later  applied  in  Section  3  to  the  SSM/T-1,2. 

2.1  General  OSSE  strategy 

There  are  four  components  common  to  any  OSSE: 

1)  A  four  dimensional  reference  atmosphere,  often  called  the  nature 
run.  This  is  considered  to  be  the  "TRUTH". 

2)  A  sampling  procedure  to  obtain  observations. 

3)  A  data  assimilation  system,  composed  of  a  forecast  model  and 
analysis  procedure. 

A)  A  quantitative  verification  procedure. 

Usually,  the  nature  run  is  simply  a  long  forecast  made  by  an  advanced  NWP 
model  or  Global  Circulation  Model.  The  more  sophisticated  the  nature  model, 
the  better.  Remotely  sensed  data  are  influenced  by  many  geophysical 
parameters,  including  sea  surface  temperature,  atmospheric  aerosol,  clouds, 
etc.  In  some  cases  these  parameters  affect  the  accuracy  of  retrievals  of 
other  parameters  or  make  such  retrievals  impossible.  These  parameters  may  be 
responsible  for  spatially  correlated  observing  errors  by  inducing  local 
geophysical  biases  in  the  retrieved  fields  (Hoffman,  1988).  These  parameters 
should  be  included  in  the  nature  run.  For  example,  SSM/T  should  provide  less 
accurate  retrieval  of  atmospheric  boundary  layer  humidity  over  the  ocean  when 
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wind  speeds  are  high  and  the  surface  emission  dominates  the  atmospheric 
emission. 

From  the  point  of  view  of  NWP,  the  most  important  characteristics  of  any 
proposed  remote  sensing  system  are  its  geographical  coverage,  horizontal  and 
vertical  resolution  and  its  error  characteristics.  In  a  simulation  study 
these  characteristics  must  be  properly  accounted  for.  These  considerations 
lead  immediately  to  a  number  of  issues  which  bear  on  the  interpretation  of  the 
results  of  OSSEs . 

For  NWP,  it  is  not  just  accuracy  of  the  measurement  which  is  important, 
the  measurement  must  be  representative  as  well.  NWP  is  really  concerned  with 
the  spatially  and  temporally  smoothed  behavior  of  the  atmosphere.  Variations 
on  the  scale  of  meters  and  seconds,  in  fact  on  the  scale  of  kilometers  and 
minutes,  are  generally  sub- grid  scale  and  are  parameterized  within  the 
model.  Consequently,  that  part  of  the  measured  signal  attributable  to  these 
scales  is  considered  to  be  noise  from  the  NWP  point  of  view.  This  source  of 
error  can  in  some  cases  be  predominant.  A  prime  example  is  radiosonde 
observations  (RAOBs).  When  two  radiosondes  are  carried  by  the  same  balloon, 
agreement  of  the  measured  quantities  is  very  good,  but  measurements  made  by 
two  radiosondes,  some  distance  apart,  do  not  agree  as  well.  One  implication 
of  this  is  that  as  models  improve  in  resolution,  this  source  of  error 
decreases.  No  existing  global  model  has  fine  enough  resolution  to  represent 
all  scales  of  motion  which  exist  in  nature.  In  fact  the  smallest  scales 
represented  by  models  are  usually  severely  damped  for  computational  reasons. 

A  method  to  unfllter  the  nature  run  was  suggested  by  Hoffman  (1988). 

The  procedures  for  simulating  data  from  the  nature  run  should  be  as 
realistic  as  possible.  The  process  of  simulating  the  observations  should  be 
sophisticated  enough  to  generate  realistic  observing  error  statistics. 

When  generating  errors  for  an  OSSE,  we  would  like  to  divide  the  process 
into  the  following  parts: 

(1)  Representational  errors.  The  nature  model  state  should  be 
unfiltered  to  restore  realistic  small  scales  to  the  model  state. 

(2)  Sampling.  The  realistic  model  state  is  interpolated  to  the 
observing  locations.  (Space  and  time  Interpolation  are  used  as  needed.) 
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(3)  Geophysical  local  bias.  Local  biases,  depending  on  the  sensor  type 
and  on  the  geophysical  parameters  of  the  realistic  model  state,  will  be 
added.  Global  biases  should  be  corrected  by  the  data  producer  and  may  be 
ignored  in  an  OSSE. 

(4)  Random  error.  The  last  error  component  is  random.  It  might  contain 
vertical  and  horizontal  correlations.  Note,  however,  that  representational 
errors  and  geophysical  local  biases  already  induce  vertical  and  especially 
horizontal  correlations.  We  feel  that  horizontal  correlations  are  mostly 
caused  by  representational  errors  and  geophysical  local  biases.  On  the  other 
hand,  vertical  errors  may  be  correlated  if  the  sensor  retrieval  algorithm 
interrelates  several  independent  observations  to  a  profile  of  retrieved 
temperature  or  other  variable. 

(5)  Sensor  filtering.  When  a  sensor  uses  a  statistical  retrieval 
method,  all  its  observations  should  be  filtered  by  projecting  onto  the 
vertical  basis  functions  which  are  used  in  the  retrieval.  This  is  also  true 
for  so-called  physical  retrieval  methods. 

Note  that  (3) -(5)  above  can  be  replaced  by  simulating  the  sensor  and  its 
retrieval  scheme  (e.g.  Atlas  et  al.,  1985).  However,  this  is  costly,  and  we 
would  prefer  instead  to  perform  selected  sensor  simulation/retrieval  st;udies 
in  order  to  define  reasonable  local  biases.  For  example,  we  might  study 
retrieval  bias  as  a  function  of  cloud  amount  for  each  retrievable  parameter. 
Such  studies  are  of  interest  in  themselves. 

Spatially  correlated  errors  are  difficult  for  an  analysis  scheme  to 
remove,  because  the  data  tend  to  corroborate  each  other.  Real  data  tend  to 
have  correlated  errors.  Even  for  radiosondes,  significant  vertical  error 
correlations  are  present.  The  data  used  In  the  present  experiments  had  data 
errors  which  are  almost  totally  uncorrelated.  As  a  result  the  errors  are  too 
easy  for  the  analysis  to  filter.  Atlas  (1988,  pers.  comm.)  reports  that  he 
undertook  some  Perfect  data  experiments  In  which  the  random  errors  were 
absent.  The  results  of  these  experiments  are  nearly  Identical  to  the  results 
of  the  experiments  which  had  random  errors. 

The  data  assimilation  system  forecast  model  and  analysis  procedures 
should  be  as  realistic  and  up  to  date  as  possible.  Impacts  of  observing 
systems  depend  to  a  certain  extent  on  the  forecast  and  analysis  methods 
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used.  In  some  cases,  there  may  be  a  severe  mismatch  between  the  analysis 
system  and  the  new  data.  For  example,  the  AFGL  analysis  system  does  not 
presently  include  a  surface  pressure  analysis  nor  does  it  use  surface  pressure 
observations  in  the  upper  air  analysis;  clearly  adding  a  new  source  of  surface 
pressure  data  would  have  no  impact.  In  general  it  must  be  noted  that 
operational  analysis  systems  have  been  tuned  for  the  data  they  normally 
receive.  The  best  use  of  a  novel  data  type  may  require  considerable  effort. 

Verification  of  OSSE  results  is  easy  because  we  have  total  knowledge  of 
the  "TRUTH" .  In  these  experiments  we  may  legitimately  use  the  word  error 
instead  of  difference  when  we  compare  an  experiment  to  the  nature  run. 
Interpretation  of  these  results  is  not  so  easy.  As  noted  above,  OSSE  results 
are  typically  too  good.  The  control  case  is  often  so  good  that  there  is 
little  room  for  positive  impact.  There  are  two  reasons  for  this:  first,  the 
forecast  model  is  more  similar  to  the  model  used  to  generate  nature  than  it  is 
to  the  real  atmosphere.  Second,  the  observational  errors  are  too  random  and 
easy  for  the  analysis  to  filter  out.  For  these  reasons  it  is  desirable  to 
calibrate  the  OSSE  results  to  OSE  results.  In  the  present  case  we  conduct  two 
OSSEs,  NOSAT  and  STATSAT,  for  which  we  have  previously  conducted  analogous 
OSEs.  We  use  only  a  very  simple  calibration  procedure  in  Section  7.- 
Basically  we  assume  OSSE  impacts  in  statistical  measures  relative  to  STATSAT 
are  proportional  to  corresponding  OSE  impacts  in  deriving  our  estimates  of 
actual  SSM/T-1,2  impacts. 

It  is  possible  to  use  a  series  of  real  analyses  for  the  reference 
atmosphere,  but  the  results  of  such  experiments  would  be  difficult  to 
interpret  for  the  following  reasons.  In  this  situation  the  "TRUTH"  is  the 
actual  atmosphere,  not  the  reference  atmosphere.  Therefore,  in  data  rich 
areas,  the  reference  atmosphere  would  agree  well  with  the  "TRUTH"  while  in 
data  voids  it  would  not.  Consequently,  simulated  observations  from  the  new 
observing  instrument  in  data  rich  areas  would  add  correct  information,  but 
have  little  impact  because  of  the  concentration  of  other  observations  already 
available,  while  simulated  observations  in  data  poor  areas  would  add  erroneous 
information,  which  would  be  carried  by  the  model  during  the  data  assimilation 
cycle  to  other  areas.  If  the  results  are  then  verified  in  data  rich  areas  we 
might  obtain  a  negative  impact  by  adding  a  new  observing  system.  Greater 
accuracy  in  the  simulated  observing  system  would  not  avoid  adding  erroneous 
data  in  data  poor  areas. 
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2.2  Sensor  simulation  methodology 

Our  general  approach  to  sensor  simulation,  which  is  specialized  to  the 
SSM/T1.2  in  Section  4  is  described  here.  Consider  a  vector  of  geophysical 
parameters.  Let  t  be  true  values  and  r  be  retrieved  values  from  a  sensor 
simulation  study.  We  will  study  the  error 

e  -  r  -  t.  (2.1) 

Our  object  is  to  model  e  in  terms  of  a  white  noise  process  so  that  we  may 
simulate  realistic  observations  using  a  pseudo-random  number  generator.  Since 
the  retrieval  has  only  a  limited  number  of  degrees  of  freedom,  we  project 
(2.1)  onto  the  empirical  orthogonal  functions  (EOFs)  of  the  retrieval  scheme, 
obtaining 


e  -  r  -  t  .  (2.2) 

Here,  for  example,  t  is  given  by 

t  -  B^(t  -  (2.3) 

where  B  is  the  matrix  of  EOFs  and  p  is  the  ensemble  average.  Our  simulated 
values  will  be  obtained  by  adding  values  of  the  modeled  e  to  values  of  t 
sampled  from  nature  and  then  transforming  back  to  physical  space. 

We  model  e  in  terms  of  geophysical  parameters,  g.  To  account  for  a 
constant  term,  we  assume  that  the  first  element  of  g  is  one.  Our  analysis 
will  yield 


e-Ag+n  (2.4) 

where  n  is  a  noise  process  with  zero  mean  and  covariance  matrix  <n  n  >. 

These  n  may  be  modeled  by 

n  -  Cc  (2.5) 
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where  C  Is  the  matrix  of  eigenvectors  of  <n  n  >  and  the  e  are  white  noise. 
Combining  (2.2),  (2.4)  and  (2.5)  yields  our  model  of  r, 

r  -  B(t  +  Ag  +  C«)  +  ft  (2.6) 

In  (2.6),  t  Is  obtained  by  projecting  the  values  from  the  nature  run  according 

to  (2.3),  g  Is  a  vector  of  geophysical  parameters  diagnosed  from  the  nature 

run  and  e  are  obtained  from  a  random  number  generator.  The  B,  A,  C  and  n 
are  constants  from  the  retrieval  algorithm  and  our  analysis  of  the  sensor 
simulation  data. 

3 .  Data  Sources 

The  data  used  In  this  study,  the  nature  run  and  slmvilated  observations 
exclusive  of  the  SSM  data  were  prepared  by  ECMWF  and  NMC  respectively  as 
described  in  a  report  by  Dey  et  al.  (1985),  which  we  summarize  below  In 
Section  3.1  and  3.2  respectively.  These  data  were  prepared  as  part  of  a  joint 
NMC- ECMWF -NASA/GLA  effort  to  conduct  WINDSAT  experiments  which  began  as  a 
result  of  a  workshop  held  at  NMC  in  February,  1983.  The  10  -  30  November 
period  was  chosen  for  this  work  because  both  ECMWF  and  GLA  had  already 
conducted  real  data  OSEs  for  that  period.  Two  TIROS  satellites  were  in 
service  during  this  period.  Further  details  of  our  OSSE  design  are  given  in 
Section  6. 

3 . 1  Nature  run 

ECMWF  generated  the  nature  run.  The  nature  run  is  simply  a  20  day. 
forecast  from  the  FGGE  Illb  analysis  produced  at  ECMWF  for  00  GMT  10  November 
1979  (Bengtsson  et  al.,  1982).  The  model  used  in  the  nature  run  forecast  was 
a  version  of  the  15  layer,  1.875*  grid  point  model  (Hollingsworth  et  al., 

1980;  Tiedtke  et  al.,  1979).  This  model  included  fairly  complete  physics  with 
a  diurnal  cycle. 

To  conserve  storage  space,  as  we  unpacked  the  gridded  nature  run  tapes  we 
interpolated  the  1.875*  grid  to  a  2.5*  grid  which  we  have  used  for  all  our 
data  sets  and  comparisons.  A  2.5*  grid  is  substantially  finer  than  the 
spectral  transform  grid  used  by  our  R30  forecast  model  and  is  therefore  more 
than  adequate  to  present  our  results.  In  fact  the  nature  run  is  rather 
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smooth,  smoother  than  many  of  our  analyses  and  forecasts,  and  the  2.5°  grid  is 
more  than  fine  enough. 

3.2  Simulated  observations 

NMC  simulated  the  FGGE  Level  lib  data  for  the  period,  in  the  NMC  fo! mat 
(Office  Note  29)  from  the  ECMWF  nature  run.  Almost. all  Level  lib  data  were 
simulated.  However  NMC  did  not  simulate  constant  level  balloon  data  (COBAL) , 
experimental  satellite  stratospheric  sounding  data  (LIMS)  and  significant 
level  data.  Later  GLA  converted  the  NMC  data  to  the  standard  FGGE  format 
(WMO,  1986).  We  received  copies  of  the  nature  run  and  FGGE  format  Level  Ilb 
data  from  GLA,  courtesy  of  R.  Atlas. 

The  simulated  standard  FGGE  Level  lib  data  were  created  by  replacing  all 
the  observed  atmospheric  variables  in  the  real  FGGE  Level  llb  data  with  values 
interpolated  from  the  nature  run  corrupted  by  adding  a  simulated  observing 
error.  Therefore  if  a  particular  radiosonde  report  is  missing  in  the  real 
data,  it  is  missing  in  the  simulated  data,  if  it  is  present  in  the  real  data, 
it  is  present  in  the  simulated  data  and  has  the  same  quality  control  marks  and 
missing  data  flags  as  the  real  observation.  This  yields  very  realistic  data 
coverage  and  quality  control  in  the  simulated  data.  However  certairt 
discrepancies  are  possible:  For  example,  CDWs  may  be  present  where  the 
moisture  field  in  the  nature  run  is  Inconsistent  with  cloudiness.  Typical 
data  coverage  by  standard  FGGE  lib  data  is  shown  in  Fig.  3.1.  The  data  points 
shown  are  those  actually  used  by  the  analysis  program,  after  a  gross  error  and 
buddy  checking  quality  control;  all  data  with  observation  times  in  a  6 -hour 
period  centered  on  the  analysis  time  are  shown.  The  radiosonde  locations  are 
those  with  height  observations  used  for  an  analysis  at  OOZ  on  25  November,  at 
the  <7—0.5  level.  All  other  plots  in  Fig.  3.1  show  data  locations  for  OOZ  on 
21  november.  Cloud  track  winds  are  most  numerous  at  low  levels  (a  -  .86  and 
.72)  and  upper  levels  (a  —  .27  and  .22),  whereas  aircraft  wind  reports  are 
concentrated  at  upper  levels  (a  -  .27  and  .22).  Coverage  by  the  TOVS  data, 
shown  here  for  height  observations  at  the  a  -  0.5  level,  illustrates  the 
satellite  tracks  and  the  fact  that  retrievals  over  land  were  not  used. 

The  value  of  the  nature  run  at  an  observing  location  is  determined  by 
spatially  interpolating  the  nature  run  at  the  closest  synoptic  time  (00,  06, 

12  or  18  GMT).  The  vertical  interpolation  is  linear  in  ln(p)  and  the 
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horizontal  Interpolation  is  quadratic  in  latitude  and  longitude  for  height 
(Z) ,  wind  components  (u,v),  and  temperature  (T) .  For  relative  humidity  (RH) , 
the  horizontal  Interpolations  are  linear. 

The  simulated  observational  error  which  is  added  to  the  value  of  the 
nature  run  at  the  observing  location  is  composed  of  a  bias  and  a  random 
Gaussian  error  which  is  not  correlated  with  anything  else.  The  size  of  the 
random  error,  or  observing  error  standard  deviation  (OESD)  is  appropriate  for 
the  particular  observation.  The  OESDs  depend  on  report  type,  variable  and 
pressure  level  and  are  displayed  in  Table  3.1  reproduced  from  Dey  et  al . 
(1985).  Biases  are  zero  except  for  TIROS. 

The  TIROS  biases  depend  on  retrieval  path  and  are  displayed  in  Table  3.2 
reproduced  from  Dey  et  al .  (1985)  which  is  based  on  Fig.  2  of  Schlatter 
(1981).  In  actual  practice,  the  retrieval  paths,  labeled  A,  B  and  C  are  set 
depending  on  whether  the  retrieval  was  deemed  clear,  partly  cloudy  or 
cloudy.  Note  that  the  OESDs  in  Table  3.1  for  TIROS  also  depend  on  the 
retrieval  path  (and  are  based  on  Fig.  3  of  Schlatter;  1981).  In  simulation 
the  retrieval  path  was  determined  from  the  nature  run  total  fractional  cloud 
coverage,  f,  according  to 


Path  “  •{ 


A 

B 

C 


if  f  <-  60 
if  60  <  f  <-  90 
if  90  <  f 


This  relationship  was  tuned  to  give  approximately  the  same  proportion  of  the 
different  retrieval  types  as  were  actually  observed  on  12  November  1979.  The 
nature  run  cloud  coverage  in  turn  is  deduced  from  the  nature  run  RH  as 
described  below.  Since  the  nature  run  RH  field  is  spatially  correlated,  the 
TIROS  observational  errors  will  be  also.  All  other  errors  are  uncorrelated. 

The  nature  run  cloud  fraction  at  500,  700  and  850  kPa  is  determined  using 
a  version  of  Fye's  (1978)  cloud  fraction  to  RH  conversion  algorithm.  This 
algorithm  is  tuned  to  the  ECMWF  forecast  so  that  the  nature  run  cloud 
statistics  are  reasonable.  Any  effect  of  high  cirrus  cloud  on  the  observation 
errors  is  ignored.  Layer  and  total  cloud  amounts  are  then  calculated  assuming 
random  overlap  between  individual  cloud  levels  within  the  layer.  (That  is, 
cloud  free  fractions  multiply  as  one  progresses  down  through  the  atmosphere.) 
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The  simulated  data  described  here  are  quite  complete  and  realistic,  yet 
have  two  major  failings.  First  the  observational  errors  are  uncorrelated. 

For  example,  consider  the  CDWs :  The  typical  CDW  OESD  used  is  of  the  proper 
size  (order  8  m/s) ,  but  real  CDW  errors  have  large  horizontal  correlations  due 
to  height  assignment  errors  which  are  responsible  for  the  largest  part  of  the 
OESD.  In  reality,  then,  the  CDW  errors  are  not  much  reduced  by  the  filtering 
of  the  analysis  procedure.  In  the  simulation  experiments  described  here,  on 
the  other  hand,  the  analysis  is  able  to  average  out  the  CDW  errors  very 
effectively  because  they  are  uncorrelated  and  the  observations  are  dense. 
Second  the  nature  run  has  little  energy  in  the  smallest  scales.  Small  scale- 
energy  present  in  the  real  atmosphere,  must  be  considered  part  of  the 
observational  error  and  thereby  induces  spatially  correlated  errors  for  all 
observations.  The  absence  of  this  source  of  error  also  contributes  to  the 
(unrealistic)  ease  with  which  the  analysis  averages  out  the  observational 
errors  in  the  simulation  experiments. 
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Table  3.1  OESDs  Used  to  Simulate  FGGE  Ilb  Data  for  November 
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SFCB  -  Surface  buoy  COB  -  COBOL  •  P  -  pressure  (in  mb) 

RAD  -  Radiosonde  SAT  -  Satellite  cloud  tracked  winds  TIRAB  -  TIROS  A&B  retrievals 

DRP  -  Dropsonde  WSAT  -  WINDSAT  TIRC  -  TIROS  C  retrievals 

PIB  -  PIBAL 


Table  3.2  Biases  Used  to  Simulate  TIROS  Temperature  Retrievals  for  November 


Retrieval  Method 
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13 


LflT 


Figure  31  Data  coverage  by  standard  FGGE  Ilb  obset^vat  ons .  Raob  heights  at 
a-O.b  (a),  cloud  track  winds  at  a-0.86  and  0.72  (b) ,  and  o»0.27 
and  0.22  (c),  aircraft  winds  at  a-0.27  and  0.22  (d) ,  and  TOVS 
heights  at  a-0 .5  (e) . 
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Figure  3.1  (continued). 


U.  Sensor  Simulation  Study  of  SSM/T  D-Matrlx  Retrievals 

Application  of  the  methodology  described  In  Section  2.2  above  to  simulate 
the  T-1  and  T-2  retrievals  required  several  processing  steps.  These  steps 
which  include  manipulation  and  extraction  of  the  nature  run  data,  radiative 
transfer  simulations  based  on  the  extracted  nature  run  data,  and  statistical 
analysis  of  the  simulation  data  are  described  in  greater  detail  below. 

A.l  Generating  dependent  and  independent  data  sets 

The  first  step  involved  generation  of  two  data  sets  containing  selected 
temperature  and  relative  humidity  profiles  culled  from  subsets  of  the  nature 
run.  The  temperature  and  humidity  data  sets,  labeled  (Dependent)  Set  1  and 
(Independent)  Set  2,  were  the  basis  for  conducting  statistical  retrievals  and 
generating  retrieval  error  statistics  which  were  used  to  simulate  retrievals 
for  the  OSSEs. 

A  particular  objective  of  our  work  is  to  realistically  reproduce  the 
error  characteristics  of  the  millimeter/microwave  observing  system,  including 
its  spatial  attributes.  Some  type  of  retrieval  classification  is  necessary 
because  geophysical  parameters  such  as  cloudiness  or  surface  type  directly 
affect  the  quality  of  the  retrieval.  The  categories  used  to  construct' the 
matrix  for  the  statistical  retrievals  (the  so  called  D-Matrix)  from  Set  1  and 
retrieval  error  statistics  from  Set  2  are  indicated  in  Table  4.1.  Note  that 
the  number  of  categories  for  water  vapor  is  twice  that  for  temperature  since 
cloudiness  has  little  effect  on  temperature  retrievals.  We  anticipated  that 
having  more  retrieval  verification  categories  would  induce  horizontal 
correlations  in  the  OSSE  retrievals. 

Sampling  of  the  nature  run  was  conducted  so  that  each  category  contained 
approximately  200  to  700  profiles.  Set  1  profiles  were  extracted  from  nature 
run  output  valid  at  00  GMT  on  17  November  and  at  12  GMT  on  26  November. 
Statistical  retrievals  were  performed  within  the  actual  time  period  chosen  for 
the  assimilation  experiments  at  12  GMT  on  21  November  and  0  GMT  on  22 
November.  Throughout  the  experiment  only  oceanic  profiles  were  used  since  the 
impact  of  the  data  would  be  greatest  in  these  regions. 

In  order  to  account  for  the  effects  of  cloud  and  ocean  surface  roughness, 
40%  of  all  profiles  were  randomly  flagged  as  cloudy  and  30%  of  all  profile.s 
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were  flagged  as  having  high  surface  wind.  These  percentages  are  consistent 
with  global  wind  (SEIASAT)  and  cloud  statistics,  (personal  communication, 

R.  Atlas,  1988;  Newell  et  al . ,  1974)  and  no  attempt  was  made  to  assign  cloud 
and  wind  flags  based  on  actual  geophysical  features.  During  retrieval 
simulations,  however,  surface  wind  and  relative  humidity,  acting  as  a 
surrogate  for  cloud  cover,  were  used  for  classification.  To  insure  realistic 
representativeness  of  the  statistics  all  conditions  are  sampled  rather  than  a 
random  distribution. 

4.2  Preprocessing  geophysical  profiles 

The  nature  run  temperature  data  are  on  12  mandatory  pressure  levels 
ranging  from  50  mb  down  to  1000  mb,  while  the  relative  humidity  data  is 
specified  on  the  lowest  6  mandatory  levels  (1000  -  300  mb) .  Since  the 
numerical  model  we  have  used  to  simulate  radiative  transfer  in  the  milli- 
meter/microwave  spectral  region,  RADTRAN  (Falcone  et  al . ,  1982),  requires 
temperature  and  water  vapor  information  on  64  pressure  levels  from  1  mb  to 
1000  mb,  preprocessing  of  the  extracted  profiles  was  performed.  Temperature 
profiles  were  interpolated  vertically  (linear  in  In  p)  between  1000  mb  and 
50  mb.  Above  50  mb  climatology  from  the  U.S.  Standard  Atmosphere  Supplements 
(1966)  was  used.  Relative  humidity  profiles  were  converted  to  water  vapor 
density  and  then  interpolated  vertically  between  1000  mb  and  300  mb,  with 
climatology  from  the  Phillips  radiosonde  set  (Phillips  et  al.,  1988)  used 
above  300  mb. 

4.3  Computing  brightness  temperatures  from  geophysical  profiles 

Radiative  transfer  simulations  were  run  for  each  of  the  12  instrument 
channels  indicated  in  Table  4.2  (Falcone  and  Isaacs,  1987).  The  brightness 
temperature  at  frequency  v,  T^(*') ,  was  evaluated  from  (Isaacs  et  al.,  1988): 


Tb(.)  - 


€  T  +  (1-6  )  (T  r;(0)  ■'T(p)dr;j 


s  s 


C  1/ 


Ps 
.  0 


T(p)dr^  (4,1) 


where 


r^(p)  -  exp  [- 


P 

k(i/,p')dp'/M] 

0 


(4.2) 
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and 


r^(p)  -  exp  [- 


®  k(i/.p')dp'//il 
P 


(^.3) 


Here  n  is  the  cosine  of  the  path  zenith  angle  and  are  the  upward  and 
downward  transmission  functions,  and  is  the  cosmic  background  tempera¬ 
ture.  Since  all  simulations  were  performed  for  nadir  fields-of-view,  the 
zenith  angle  cosine,  in  (4.2)  and  (4.3)  above  was  set  equal  to  unity. 

The  nonscattering  radiative  transfer  equation  (4.1)  Is  used  for  all 
simulations  since  scattering  by  cloud  size  droplets  is  negligible.  The  cloud 
model  drop  size  distributions  (Table  4.3)  have  mode  radii  of  less  than  10  fim 
with  appreciable  numbers  of  drops  falling  off  rapidly  beyond  100  fim.  The 
single  scattering  albedo  for  such  drops  at  183  GHz  is  about  10'^.  Generally, 
a  single  scattering  albedo  of  0.05  or  greater  is  required  before  a  significant 
scattering  source  function  is  attained.  At  183  GHz,  this  implies  significant 
numbers  of  drop  sizes  of  a  few  millimeters  or  more,  i.e.  precipitation  sized 
drops.  Due  to  the  nature  of  the  cloud  drop  size  distributions  treated,  there 
are  few  such  large  droplets  (OeirmendJ ian,  1975). 

Profiles  of  atmospheric  absorption  coefficients,  k(i/,p),  for  each 
frequency  of  interest  were  calculated  using  the  RADTRAN  simulation  algorithm, 
and  these  provided  the  transmission  profiles  (4.2),  (4.3)  necessary  to 
evaluate  channel  brightness  using  (4.1).  The  channel  set  and  design  noise 
equivalent  brightness  temperatures  (NEAT' s)  for  the  SSM/T-l  and  T-2  sensors 
are  given  in  Table  4.2.  For  Set  1  profiles  no  noise  was  added  to  computed 
brightness  temperatures  since  these  were  only  used  to  construct  D-Matrices. 
During  simulations  using  Set  2  profiles  Gaussian  noise  with  a  mean  of  zero  and 
standard  deviation  equal  to  the  channel  NEAT  was  added  to  each  brightness 
temperature . 

The  procedure  described  above  was  applied  to  clear  sky  simulations.  To 
simulate  the  effect  of  cloud  within  the  fleld-of-vlew  of  the  radiometer,  a 
suitable  set  of  cloud  models  was  incorporated  within  the  evaluation  of  the 
atmospheric  attenuation  profiles.  The  cloud  models  chosen  were  taken  from  the 
AFGL  FASCODE  model  (Falcone  et  al . ,  1979).  Mode  radii,  cloud  liquid  water 
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content  (LWC) ,  and  the  vertical  extent  of  the  five  cloud  models  are  given  in 
Table  A. 3. 

Cloud  attenuation,  a(i/,T,z),  was  calculated  using  the  Rayleigh 
approximation  (van  de  Hulst,  1957): 

2 

a(i/,T,z)  -  6x1/  Im  (-  )  LWC  (4.4) 

2 

where  »/  is  the  frequency  in  wavenumbers,  m  is  the  complex  index  of  refraction 
of  water  evaluated  using  Ray's  (1972)  empirical  expression,  and  LWC  is  the 
cloud  liquid  water  content.  Falcone  et  al. ,  (1979)  have  shown  that  the  use  of 
the  Rayleigh  approximation  Instead  of  the  full  Mie  theory  for  these  cloud 
models  is  justified  for  frequencies  less  than  about  300  GHz. 

Clouds  were  incorporated  Into  the  brightness  temperature  simulation 
process  by  randomly  selecting  from  among  one  of  the  cloud  types  with 
associated  LWC  given  in  Table  4.3  or  choosing  the  option  that  the  simulated 
scene  was  clear.  Cloud  attenuation  calculated  from  (4.4)  was  then  added  to 
the  clear  sky  absorption  for  each  atmospheric  layer  within  the  vertical  domain 
specified  for  the  selected  cloud  type  model.  The  effect  of  partially  filled 
fields-of -view  was  then  treated  by  multiplying  the  local  cloud  attenuation  by 
a  random  number  between  0  and  1.  It  was  assumed  that  all  cloud  types  were 
equally  probable,  and  that  only  a  single  cloud  layer  was  present.  The  second 
assumption  tends  to  underestimate  cloud  Impact,  but  we  note  that  these 
assumptions  could  be  relaxed  given  appropriate  cloud  coverage  statistics  and 
frequency  of  occurrence  data.  Given  that  beam  filling  is  particularly 
unlikely  for  cumuliform  clouds  with  horizontal  extents  of  a  few  kilometers 
this  treatment  of  the  f ield-of-view  should  be  more  accurate.  Since  the 
profiles  of  temperature  and  water  vapor  were  given  at  constant  pressure  levels 
rather  than  at  constant  heights  as  are  the  cloud  models,  cloud  vertical  extent 
had  to  be  interpolated  to  the  appropriate  pressure  level  within  individual 
profiles.  Within  cloudy  layers  relative  humidity  was  not  adjusted  to 
saturation. 

Finally,  for  millimeter  and  microwave  radiation  surface  emissivity  is 
typically  less  than  1.0  and  retrieval  results  are  highly  sensitive  to 
characteristics  of  the  underlying  surface  (Isaacs  and  Deblonde,  1987).  During 
simulations  we  allowed  the  emissivity  to  vary  with  sensor  and  surface  type 
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according  to  Table  4.4.  These  are  consistent  with  the  RADTRAN  surface 
emissivity  models  given  by  Isaacs  et  al.  (1989). 

4.4  D-Matrix  and  statistical  retrievals 

A  variety  of  possible  retrieval  approaches  could  be  applied  to  the 
183  GHz  water  vapor  retrieval  problem.  These  include  both  physically  based 
methods  and  those  relying  on  statistical  principles.  These  are  summarized  in 
Table  4.5  (see  Isaacs,  1987  for  further  details).  For  this  study  we  have 
adopted  a  simple  statistical  approach. 

The  specific  statistical  inversion  or  D-Matrlx  approach  to  atmospheric 
parameter  retrievals  has  been  described  by  various  investigators  (Rodgers, 
1976;  Smith  and  Woolf,  1976;  Gaut  et  al.,  1975;  Isaacs  et  al . ,  1985).  It  is 
essentially  a  least  squares  approach  which  finds  the  most  likely  combination 
of  atmospheric  parameters  which  yields  the  set  of  observed  radiometric  data. 
The  method  used  here  utilizes  an  eigenanalysis  of  the  observed  data  and 
retrievable  parameter  covariance  matrices.  Briefly,  an  Individual  retrieval 
of  temperature  or  relative  humidity  on  k  vertical  levels  from  n  channels  of 
brightness  temperature  measurements  is  represented  by  the  matrix  equation: 

p  -  Db 

with 

D  -  (Pb'^)  ($A''«>'^) 

where 

p  is  a  vector  containing  the  retrieved  profile 

b  is  a  vector  containing  n  brightness  temperatures 

is  the  retrievable  atmospheric  profile  at  k  levels  for  s 
samples 

Bj^g  are  brightness  temperatures  for  n  channels  and  s  samples 

T 

are  eigenvectors  of  BB 

A  is  a  diagonal  matrix  whose  elements  are  the  corresponding 
eigenvalues 

By  setting  some  of  the  diagonal  elements  of  A  to  zero,  the  method  allows  one 
to  retain  only  those  EOFs  whose  corresponding  eigenvalues  are  relatively 
large.  This  yields  a  more  stable  matrix  inversion  and  reduces  the  effect  of 
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noise  in  the  brightness  temperature  data.  In  this  study,  then,  the  tempera¬ 
ture  and  relative  humidity  profiles  and  associated  brightness  temperatures 
from  Set  1  were  combined  to  construct  D-Matrlces  for  each  of  the  categories 
shown  in  Table  A.l.  An  eigenvector  was  retained  only  if  the  ratio  of  its  cor¬ 
responding  eigenvalue  to  the  largest  eigenvalue  was  greater  than  1.0  x  10"^. 

Once  the  necessary  D-Matrices  were  created,  error  statistics  were 
calculated  by  simply  retrieving  relative  hvunldity  and  temperature  from  Set  2 
radiances  using  the  Set  1  D-Matrices  and  comparing  the  retrieved  profiles  with 
the  true  profiles.  Retrievals  of  both  variables  were  done  separately  in  that 
humidity  was  first  retrieved  using  both  SSM/T-1  and  SSM/T-2  data,  followed  by 
temperature  retrievals  using  the  same  radiometric  data. 

For  each  error  statistic  category  in  Set  2  the  retrieval  or  observing 
error  standard  deviation  (OESD)  at  each  level,  the  mean  retrieval  error 
(bias) ,  and  the  interlevel  correlations  of  the  retrieval  error  were 
computed.  Subsequently,  EOFs  of  the  vertical  correlations  were  also  computed; 
these  were  needed  to  reconstruct  the  appropriate  vertical  errors  during  the 
retrieval  simulations.  A  more  detailed  discussion  of  the  actual  retrieval 
errors  which  were  used  during  simulation  is  contained  in  Section  5. 

4.5  Simulation  of  D-Matrix  retrievals 

To  conduct  the  OSSE  the  following  sequence  of  steps  was  repeated  for  the 
location  of  each  simulated  retrieval: 

•  Compute  the  next  satellite  field-of-view  (FOV)  location  and  time, 

•  Interpolate  the  nature  run  data  to  the  FOV  location. 

•  Based  on  predefined  criteria,  determine  the  geophysical  category  of 
profile . 

•  Perturb  the  Interpolated  profile  based  on  OESD,  bias  and  vertical 
error  correlations  for  that  category. 

Simulations  were  created  for  the  entire  assimilation  period  (6  Z 
18  November  1979  -  0  Z  25  November  1979).  Retrievals  were  organized  in  6  hour 


22 


intervals,  centered  on  the  synoptic  times  to  accommodate  the  AFGL  statistical 
analysis  system  (Norquist,  1986). 

Computation  of  the  satellite  surface  FOV  requires  specification  of 
(1)  the  orbital  characteristics  of  the  sensor-bearing  platform,  and  (2)  the 
scanning  attributes  of  the  sensor  itself.  The  orbital  parameters  of  the  DMSP 
polar  orbiting  satellite  are  well  documented  and  we. have  used  these  in  our 
simulations.  Although  the  operational  T-1  and  proposed  T-2  instruments  have 
different  scan  patterns  we  have  assumed  for  consistency  that  both  scan 
patterns  are  identical  so  that  measurements  from  one  sensor  are  coincident 
with  another.  This  eliminates  an  additional  preprocessing  ambiguity  in  which 
smaller  T-2  footprints  would  have  to  be  combined  and  interpolated  to  the 
larger  FOV  of  the  T-1  sensor.  Orbital  and  scan  parameters  as  well  as  details 
of  the  FOV  computation  used  in  our  simulations  are  described  in  Appendix  A. 
Note  that  although  scanning  characteristics  associated  with  the  T-2  sensor  are 
included  in  Appendix  A  they  were  not  used  during  the  OSSE. 

Once  the  location  of  the  current  sounding  location  is  computed,  it  is 
verified  against  a  data  base  of  surface  type  and  any  FOV  falling  over  land  is 
rejected.  For  an  accepted  FOV,  nature  run  data  is  interpolated  bilinearly  in 
space  and  linearly  in  time  from  the  nearest  two  nature  run  time  levels. to 
create  the  unperturbed  or  "true"  profile. 

The  resulting  data  pattern  for  a  typical  6-hour  period  (in  this  case  OOZ 
on  25  November)  is  shown  in  Fig.  4.1. 

Assignment  of  the  profile  into  one  of  several  geophysical  classes  is 
based  on  a  number  of  criteria: 

Surface  type  (ice/ocean):  From  an  NMC  data  base  which  reflects  November 

climatology,  each  gridpoint  on  the  2.5*  x  2.5*  grid  is  flagged  as  either 
land,  ice  covered,  or  ocean. 

Geographic  location:  Points  poleward  of  approximately  30°  which  are  not  ice 
covered  are  flagged  as  midlatitude  profiles.  All  other  profiles  are 
considered  tropical.  To  allow  for  a  smooth  transition  from  tropical  to 
midlatitude  statistics,  the  crossover  point  for  each  satellite  orbit  was 
taken  to  be  30°  plus  a  random  Gaussian  perturbation. 
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Sea  surface  roughness:  Assuming  that  the  1000  mb  wind  in  the  unperturbed 
profile  is  equal  to  the  surface  wind,  profiles  with  a  horizontal  wind 
speed  less  than  or  equal  to  7  ms*^  are  assigned  a  low  wind/calm  ocean 
flag.  Those  with  a  higher  surface  wind  are  assigned  the  high  wind/rough 
ocean  flag. 

Cloud  cover:  Methods  to  infer  cloud  cover  from  relative  humidity  data  (and 
vice-versa)  vary  in  complexity.  The  approach  we  have  adopted  utilizes 
the  Tibaldi  scheme  (Norquist,  1988)  which  is  actually  a  means  of  convert¬ 
ing  cloud  cover  data  to  relative  humidity  in  four  atmospheric  layers.  We 
have  adapted  the  method  to  solve  the  inverse  problem,  and  since  the 
scheme  then  yields  a  percent  cloud  cover  in  each  vertical  layer  we 
convert  this  to  a  binary  cloud/no  cloud  assignment.  If  the  number  of 
relative  humidity  levels  in  which  the  computed  cloud  cover  exceeds  50%  is 
greater  than  3,  then  the  profile  is  cloudy.  In  all  other  cases  the 
profile  is  clear. 

Finally,  the  profile  derived  from  nature  run  output  was  perturbed  by  an 
amount  which  depends  upon  the  error  statistics  associated  with  the  geophysical 
category.  Given  an  unperturbed  profile  of  atmospheric  parameters,  P,  Che 
perturbation  is  of  the  form: 

P  -  P  +  E 

where 

P  is  the  unperturbed  profile  obtained  from  the  nature  run 

E  is  the  retrieval  error 

P  is  the  retrieved  profile. 

The  total  retrieval  error  at  any  vertical  level  k  may  be  written  as 

^  k 

where 

is  the  systematic  component 

# 

E  is  a  random  component. 
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For  each  profile  to  be  retrieved  the  systematic  component  is  the  same  for  each 
category  and  level.  The  random  component  is  computed  from: 


\  -  r®lk 


with 


i-1 


-  ^k  ‘^k 


where 


is  the  value  of  the  EOF  at  level  k  determined  from 

eigenanalysis  of  the  vertical  error  correlations 


is  the  retrieval  error  standard  deviation 

e  is  a  random  number  from  a  Gaussian  distribution  with  a  mean  of 
0  and  a  standard  deviation  of  1 


N  is  the  number  of  EOFs. 

Since  the  error  correlation  EOFs,  the  error  standard  deviation,  and  the 
systematic  error  are  different  for  each  geophysical  category  horizontal  and 
vertical  retrieval  error  statistics  should  be  reproduced  with  greater 
fidelity . 
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Table  4.1  Geophysical  categories  used  to  generate  retrieval  error  statistics.  Also 
shown  are  the  number  of  samples  in  each  category. 


Parameter 

Set  1 
Training 

Sub -Sample 

Number  of 
Profiles 

Set  2 

Geophysical  Independent  Number  of 

Category  Sub- Sample (s)  Profiles 

T 

Ice 

438 

1 

Ice 

438 

Ocean 

536 

2 

Calm  ocean.  Tropics 

459 

3 

Calm  ocean,  Mid-Latitudes 

493 

4 

Rough  ocean.  Tropics 

583 

5 

Rough  ocean,  Mid-Latitudes 

629 

RH 

Ice,  Clear 

430 

1 

Ice,  Clear 

437 

Ice,  Cloudy 

452 

2 

Ice,  Cloudy 

440 

Ocean,  Clear 

430 

3 

Calm  ocean.  Tropics,  Clear 

425 

4 

Calm  ocean,  Mid-Latitudes,  Clear 

523 

5 

Rough  ocean.  Tropics,  Clear 

435 

6 

Rough  ocean,  Mid-Latitudes,  Clear 

293 

Ocean,  Cloudy 

427 

7 

Calm  ocean.  Tropics,  Cloudy 

351 

8 

Calm  ocean,  Mid-Latitudes,  Cloudy 

238 

9 

Rough  ocean.  Tropics,  Cloudy 

392 

10 

Rough  ocean,  Mid-Latitudes,  Cloudy 

259 

26 


Table  4.2  Microwave  mission  sensor  characteristics 


Instrument 

Frequency 

(GHz) 

Polari¬ 
zation 
(H  or  V) 

FOV 

(km) 

Response 

NEAT 

(K) 

SSM/T-1 

50.5 

H 

200 

surface 

0.6 

53.2 

H 

200 

T  at  2  km 

0.4 

54.35 

H 

200 

T  at  6  km 

0.4 

54.9 

H 

200 

T  at  10  km 

0.4 

58.825 

V 

200 

T  at  16  km 

0.4 

59.4 

V 

200 

T  at  22  km 

0.4 

58.4 

V 

200 

T  at  30  km 

0.5 

SSM/T-2 

91.655il.25 

V 

100 

surface,  water 

vapor 

0.6 

150.00±1.25 

V 

60 

surface,  water 

vapor 

0.6 

183. 31±1 

V 

50 

water  vapor 

0.8 

183 . 31±3 

V 

50 

water  vapor 

0.6 

183.13±7 

V 

50 

water  vapor 

0.6 
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Table  4.3  Cloud  type  characteristics 


Model 

Type 

Mode  Radius 
(fua) 

Liquid  Water 
Concent  (gm"'^) 

Vertical  Extent 
(km) 

1 

Stratus 

2.7 

0.15 

0.5  -  2.0 

2 

Cumulus 

6.0 

1.00 

1.0  -  3.5 

3 

Altostratus 

4.5 

0.40 

2.5  -  3.0 

4 

Stratocumulus 

6.25 

0.55 

0.5  -  1.0 

5 

Nimbostratus 

3.0 

0.61 

0.5  -  2.5 
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Table  ■!» .  4  Surface  emissivity  used  for  various  surface  types  during  radiative 
transfer  simulations 


I  n,strunient 

Ice 

Calm  Ocean 

Rough  Ocean 

SSM/T-1 

.85 

.45 

.58 

SSM/T-2 

,95 

.  70 

.83 
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Table  4.5  183  GHz  water  vapor  profile  retrieval  techniques 


Investigators/ 

Date 

Organization  Approach 

Reference 

Gaut  et  al . , 
1975 

ERT 

Use  of  183  GHz  and 
other  channels  for 
water  vapor  and  cloud 
properties 

AFGL-TR-75-0007. 

Schaerer  and 
Wilheit,  1979 

NASA 

Profile  retrieval  simu¬ 
lation  using  5  channels 
and  iterative  method 

Radio  Science.  14. 

3,  371-375. 

Rosenkranz 
et  al. .  1982 

MIT 

Statistical  retrieval 
based  on  60  GHz  and 

183  GHz  simulations 

J.  AddI.  Meteor..  21. 
1364-1370. 

Wang  et  al . , 
1983 

NASA 

Aircraft  radiometer  data 
and  retrieval 

J .  Glim,  and  Aool . 
Meteor. .  22.  779-788. 

Kakar,  1983 

JPL 

Chahine-type  retrieval 

J.  Glim,  and  Anol. 
Meteor..  22.  1282-1289. 

Kakar  and 
Lambrigsten, 
1984 

JPL 

Statistical  correlation 
technique 

J  .  Clint .  and  Aool . 
Meteor..  23.  1110-1114. 

Isaacs  and 
Deblonde 

1985,  1987 

AER 

Statistical  retrieval 

AFGL-TR-85-0040; 

AFGL-TR-85-0095. 

Radio  Science.  22  3. 
367-377. 

Isaacs  et  al . , 
1988 

AER 

Unified  retrieval 

AFGL-TR-88-0058. 
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5.  Statistical  Analysis  of  Simulated  Sensor  Errors 

After  creating  data  sets  1  and  2  which  contained  selected  relative 
humidity  and  temperature  profiles  spanning  a  range  of  geophysical  conditions 
D-Matrix  statistics  were  computed  and  retrievals  were  conducted  as  described 
in  Section  4.  These  calculations  resulted  in  a  set  of  class  dependent 
retrieval  error  statistics  which  were  used  in  the  subsequent  observing  system 
simulations . 

5.1  Relative  humidity  retrieval  error  statistics 

Retrieval  bias 

The  systematic  error  or  bias  of  the  relative  humidity  (RH)  retrievals 
which  were  obtained  by  retrieving  set  2  profiles  with  set  1  D-Matrices  are 
shown  in  Figs.  5.1  through  5.10.  Each  vertical  profile  is  representative  of  a 
particular  geophysical  category;  the  categories  themselves  have  been  described 
in  section  4. 

Midlatitude  profiles  (Figs.  5.5,  5.6,  5.9,  5.10)  have  the  same  general 
shape  consisting  of  mostly  positive  bias  throughout  with  larger  values 
(approximately  5  to  10  percent  RH)  in  the  mid  troposphere  decreasing  to  near 
zero  at  the  surface  and  300  mb.  At  the  same  time  tropical  bias  profiles 
(Figs.  5.3,  5.4,  5.7,  5.8)  tend  to  have  a  structure  with  negative  values  in 
the  lower  troposphere  changing  to  positive  values  at  higher  levels.  Profiles 
for  polar  ice  surface  retrievals  show  little,  if  any,  bias  at  the  six 
retrieval  levels.  And  in  all  cases  retrieval  bias  is  generally  small,  in  the 
range  of  5  to  10  percent,  but  occasionally  increasing  to  near  20  percent  (see 
Class  3,  4,  and  7  profiles). 

Finally,  there  appears  to  be  only  a  slight  difference  between 
corresponding  clear  and  cloudy  profiles,  the  only  noticeable  difference  being 
a  slight  negative  shift  (about  5  percent)  of  the  cloudy  profiles  to  more 
negative  (or  less  positive)  values.  This  may  be  a  function  of  the  cloud 
models  used  and  f ield-of -view  beam  filling  effects,  however  the  fact  that 
separate  D-Matrices  for  clear  and  cloudy  conditions  were  used  ought  to 
minimize  the  cloud  signature  in  the  bias  profiles. 
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Interpreting  the  nature  of  these  systematic  errors  is  not  entirely 
straightforward.  The  relationship  of  variations  in  observed  brightness 
temperatures  to  variations  in  water  vapor  amounts  under  conditions  of  non-unit 
surface  emissivity  is  highly  dependent  upon  the  value  of  the  emissivity  and 
the  structure  of  the  low  level  temperature  profile,  as  well  as  the  water  vapor 
amount  itself.  Moreover,  the  water  vapor  retrieval  problem  is  inherently  non¬ 
linear  because  variations  in  w.^ter  vapor  concentrations  directly  affect  the 
atmospheric  transmission  and  therefore,  the  vertical  weighting  functions. 

This  implies  that,  for  example,  an  increase  in  water  vapor  in  the  1000  to  850 
mb  layer  may  produce  either  a  positive  or  a  negative  change  in  the  T-2  channel 
radiances  and  there  is  no  clear  way  to  state  a  priori  what  the  sign  of  this 
change  will  be. 

Retrieval  error  standard  deviation 

The  OESD  profiles  for  relative  humidity  retrievals  are  shown  in 
Figs.  5.11  through  5.20.  These  profiles  generally  contain  larger  errors  at 
700  or  500  mb  and  lower  errors  at  the  surface.  This  is  at  variance  with 
typical  statistical  retrieval  errors  which  tend  to  have  largest  errors  in  near 
surface  layers. 

The  range  of  errors  in  the  OESD  curves  is  considerably  greater  than  the 
corresponding  bias  profiles  with  values  between  8  and  25  percent.  This  means 
that,  on  average,  random  errors  will  tend  to  dominate  the  systematic  error 
component  and  any  error  correlations  in  the  horizontal  between  simulated 
retrievals  will  tend  to  be  small.  This  was  indeed  the  case  when  resulting 
retrieval  error  correlations  of  the  OSSE  were  computed  as  a  function  of 
distance . 

Geographical  dependence  of  the  OESD  profiles  is  also  evident.  The 
essential  feature  here  is  that  extratropical  profiles  (Figs.  5.15,  5.16,  5.19, 
5.20)  have  larger  values  than  corresponding  tropical  profiles  (Figs.  5.13, 

5.14,  5.17,  5.18).  This  is  most  likely  a  reflection  of  the  greater 
atmospheric  variance  present  in  the  midlatitudes  than  in  the  tropics.  Since 
the  retrieval  scheme  operates  in  a  least  squares  sense,  it  tends,  in  the  mean, 
to  retrieve  profiles  reflective  of  the  average  atmospheric  conditions  which 
were  used  to  construct  the  particular  D-Matrix.  Thus  statistical  retrievals 
of  profiles  from  a  population  with  higher  variance  will  result  in  higher  rms 
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error  values.  This  tropical/extratroplcal  distinction  has  probably  been 
blurred  somewhat  by  binning  profiles  from  the  northern  hemisphere  extratropics 
together  with  southern  hemisphere  extratropical  profiles  because  conditions  in 
the  November  southern  hemisphere  (emerging  summer)  are  likely  to  be  more 
similar  to  tropical  conditions  than  those  in  the  northern  hemisphere 
atmosphere  (emerging  winter) . 

Finally,  the  effect  of  cloudiness  is  quite  evident  in  these  profiles.  A 
comparison  of  clear/cloudy  category  pairs  shows  that  cloudy  profiles  have 
errors  2  to  4  percent  greater  than  their  clear  counterparts.  Moreover,  the 
effect  is  greatest  at  levels  below  500  mb  where  most  of  the  clouds  in  the 
forward  model  are  found. 

Vertical  error  correlations 

Vertical  correlations  of  the  retrieval  errors  were  computed  for  all 
categories  and  two  examples  are  shown  here.  Fig  5.21  shows  the  correlation  of 
errors  at  1000  mb  with  all  other  levels  for  Class  1  retrievals.  In  this 
example  the  errors  at  levels  300  mb  through  850  mb  are  only  weakly  correlated 
with  that  at  the  surface.  Fig.  5.22  shows  the  interlevel  correlations  with 
respect  to  the  500  mb  level.  This  is  an  example  of  a  more  typical  bimodal 
structure  with  relatively  high  positive  correlations  at  adjacent  levels  and 
smaller  anticorrelations  at  more  distant  levels. 

5.2  Temperature  retrieval  error  statistics 

Retrieval  bias 

Profiles  of  temperature  retrieval  systematic  errors  for  five  categories 
are  shown  in  Figs.  5.23  through  5.27.  The  most  significant  differences  among 
the  profiles  depend  on  latitude.  For  example,  both  midlatitude  profiles 
(Figs.  5.25  and  5.27)  are  characterized  by  negative  biases  of  order  -0.5 
to  -1.0  K  below  300  mb  and  smaller  positive  biases  above  this  level.  On  the 
other  hand,  profiles  obtained  from  tropical  retrievals  (Figs.  5.24  and  5.26) 
have  generally  positive  biases  of  order  0.5  to  1.0  K  in  the  lower  to  middle 
troposphere  with  negative  errors  at  higher  levels.  These  bias  values  in  the 
lower  atmosphere  are  somewhat  puzzling.  Tlie  statistics  used  to  construct  the 
retrieval  matrix  were  an  amalgam  of  both  tropical  (relatively  warm)  and 
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midlatitude  (relatively  cold)  profiles.  Since  the  single  D-Matrlx  was  used 
one  might  expect  that  there  would  be  a  tendency  for  the  retrieved  lower 
troposphere  temperatures  to  be  too  high  in  midlatitude  cases  and  too  cold  in 
tropical  cases.  The  results  indicated  here  may  again  be  a  reflection  of  the 
fact  that  northern  hemisphere  winter  and  southern  hemisphere  summer  profiles 
were  mixed  in  computing  the  midlatitude  statistics.  As  is  the  case  with  RH 
retrievals,  the  effect  of  varying  sea  surface  roughness  (surface  emissivity) 
upon  the  error  profiles  is  juite  weak. 

Retrieval  error  standard  deviation 

All  the  temperature  OESD  profiles  contain  large  values  near  the  surface 
and  also  near  the  tropopause.  This  agrees  with  results  found  in  previous 
efforts  to  retrieve  temperature  using  microwave  channels  (Isaacs  et  al .  , 
1988).  And  the  latitude  dependence  of  the  vertical  error  structure  is  again 
consistent  with  notions  of  atmospheric  variability  which  applied  to  the  RH 
retrievals.  That  is,  errors  are  higher  by  approximately  1.0  K  in  the 
midlatitude  categories  (Figs.  5.30  and  5.32)  than  in  the  tropical  categories 
(Figs.  5.29  and  5.31).  For  all  profiles  the  range  of  values  is  from  1.0  to 
3.0  K  at  all  levels. 

Vertical  error  correlations 

Lastly,  the  interlevel  correlations  of  the  temperature  retrieval  error 
were  computed  and  several  examples  of  correlations  for  Class  1  are  shown  in 
Figs.  5.33  through  5,35.  Each  curve  has  a  clearly  defined  structure; 
retrieval  errors  are  positively  correlated  at  nearby  levels  and  negatively 
correlated  for  levels  which  are  more  distant. 

In  sum,  these  statistics  were  compiled  by  performing  actual  D-Matrix 
retrievals  using  simulated  brightness  temperatures  corresponding  to  a  small 
subset  of  nature  run  data.  These  statistics,  in  turn,  were  an  integral  part 
of  the  OSSE  which  was  conducted  using  the  much  larger  data  set  of  the  week 
18  -  25  November  1979. 
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Figure  5.2  Relative  humidity  retrieval  bias,  Class  2. 
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Figure  5.3  Relative  humidity  retrieval  bias,  Class  3 
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gure  5.4  Relative  humidity  retrieval  bias,  Class  4 
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Figure  5.7  Relative  humidity  retrieval  bias,  Class  7. 
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Figure  5.8  Relative  humidity  retritval  bias,  Class  8 
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Figure  5.16  Relative  humidity  retrieval  error,  standard  deviation,  Class  6. 
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Figure  5.17  Relative  humidity  retrieval  error, 


22  24  26  28  30 


standard  deviation,  Class  7 
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Figure  5.18  Relative  humidity  retrieval  error,  standard  deviation,  Class  8. 


T — I — I — I — I — I — I — I — I — I — I — I — I — I  I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 


■  I  ■  1  »  I  ■  i-t'  1  ■  I  ■  1  ■  1  ■  I  ■  I  I  I  I  I  I  I  .  I  . 


53 


PRESSURE  (  MILLIBARS  ) 


D-WTRIX  R.H.  RETRIEVAL  ERROR 
CUBS  S 


Figure  5.19  Relative  humidity  retrieval  error,  standard  deviation,  Class  9 
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Figure  5.20  Relative  humidity  retrieval  error,  standard  deviation,  Class  10 
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Figure  5.23  Temperature  retrieval  bias,  Class  1. 
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Figure  5.28  Temperature  retrieval  error,  standard  deviation 


Figure  5.29  Temperature  retrieval  error,  standard  deviation,  Class  2 
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Figure  5.31  Temperature  retrieval  error,  standard  deviation, 
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Figure  5.32  Temperature  retrieval  error,  standard  deviation,  Class  5. 
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6,  OSSE  Design 

We  now  describe  the.  OSSE  we  have  conducted  which  makes  use  of  the 
sinuilated  data  described  above.  In  our  experiments  the  data  assimilation 
sYsti'in.  used  is  the  AFGL  system  which  was  used  by  Norquist  (1988)  and  Louis  et 
al  (1*88).  This  system  is  briefly  described  in  Section  6.1,  with  special 
atreiition  to  changes  made  for  the  present  experiments.  Our  verification 
strategy  is  stu'.imar i zed  in  Section  6.2. 

To  rriniroize  transient  behavior,  i.e.  spinup,  we  preceded  our  experiment.'; 
wit!,  a  spin'up  procedure  which  included  a  four  day  forecast  and  three  day  data 
as  s  i:v.  i  i  a  t  ion  as  described  in  Section  6.3.  The  experiments  themselves, 

STATSAT,  NOSAT,  SSM+TOVS  and  SSMSAT  and  the  companion  real  data  experiments 
are  described  in  Section  6.4.  The  real  data  OSEs  arc  described  in  greater 
ceta;  i  by  I.oui  s  et  al  .  (1988)  . 

6  1  .AFGL  forecast  and  analysis  system 

Each  simulated  data  assimilation  experiment  described  here  consists  of 
one  assimilation  run  for  seven  days  and  three  forecasts,  each  four  days  in 
Length.  The  real  data  assimilation  experiments  each  consist  of  two  of  these 
assimilation  runs.  Each  assimilation  run  consists  of  a  series  of  assimilation 
cycles,  and  each  cycle  in  turn  is  made  up  of  a  6 -hour  forecast  that  serves  as 
a  first  guess  for  the  analysis,  an  optimum  interpolation  analysis  which 
combines  the  first  guess  fields  with  the  observations,  and  a  nonlinear  normal 
mode  initialization  of  the  analysis.  The  initialized  analysis  is  the  starting 
point  for  the  next  6-hour  forecast,  which  is  then  used  as  the  first  guess  of 
the  subsequent  assimilation  cycle.  The  forecast  model  used  for  the  6-hour 
foreca.st  is  a  complete  global  spectral  model  (GSM).  This  model  is  also  used 
to  produce  forecasts  out  to  4  days  starting  from  days  3,  5,  and  7  of  the 
assimilation  runs. 

6.1.1  Analysis 

The  AFGL  Statistical  Analysis  Program  (ASAP)  (Norquist,  1986,  1988)  was 
de/eloped  from  the  NMC  multivariate  optimal  interpolation  (01)  procedure  as 
d.  scribed  by  Bergman  (1979)  and  by  McPherson  et  al.  (1979).  The  ASAP  01  is  a 
miltivariate  analysis  of  height  and  wlrd  components  and  a  univariate  analy.si.s 
of  relative  humidity,  both  in  model  sigma  layers.  The  correc t ion.<='  for  an 
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analysis  grid  point  are  weighted  sums  of  surrounding  observation-minus- first 
guess  residuals.  The  equations  for  these  weights  as  well  as  the  computation 
of  the  horizontal  and  vertical  correlation  functions  follow  Bergman  (1979). 

The  analysis  error  evolves  according  to  simple  rules  (Norquist,  1986).  The 
great  circle  distance  method  for  correlation  functions  equatorward  of  70°C 
latitude  is  included  as  described  by  Dey  and  Morone  (1985)  without  changing 
the  Bergman  formulation  (including  map  factor)  for  latitudes  poleward  of  70° 
latitude.  The  analysis  takes  place  in  the  sigma  coordinates  of  the  model  on  a 
Gaussian  grid  of  62  x  61  latitude-longitude  points. 

Data  used  by  the  height-wind  analysis  include  Type  1  observations  (radio¬ 
sondes,  pibals,  etc.),  Type  2  observations  (aircraft).  Type  4  observations 
(satellite  retrieved  temperatures  or  thicknesses)  and  Type  6  observations 
(cloud  drift  winds  (CDWs)).  The  Type  3  surface  observations  are  not  used  at 
all.  This  implies  that  satellite  "heights"  are  anchored  only  by  the  6  h  fore¬ 
cast  in  regions  where  radiosondes  are  absent.  In  all  experiments,  except  for 
the  three-day  preliminary  assimilations,  the  CDW  data  were  combined  (i.e. 
locally  averaged)  into  "super-obs".  There  are  two  principal  reasons  for  doing 
this:  First,  to  limit  the  total  number  of  observations,  so  that  computer 

memory  restrictions  are  not  exceeded,  and  second,  the  CDW  errors  are.  strongly 
correlated  horizontally  because  the  main  error  is  due  to  height  assignment. 
Satellite  temperature  profiles  are  not  used  over  land  in  any  of  the 
experiments . 

Since  the  statistical  models  used  by  the  01  are  never  exact  in  practice, 
we  decided  for  convenience  and  realism  to  leave  most  of  the  statistical  models 
in  the  01  as  they  were  for  the  real  data  OSEs  described  by  Louis  et  al . 

(1988).  The  statistical  models  and  parameters  used  are  identical  to  those 
described  by  Norquist  (1986),  which  in  turn  are  based  on  NMC  practice  as 
described  by  Dey  and  Morone  (1985).  The  01  assumes  that  radiosonde 
observational  errors  are  correla-ed  vertically  and  that  satellite  height 
observational  errors  are  correlated  vertically  and  horizontally.  A  number  of 
studies  colocating  satellite  and  radiosonde  height  data  were  performed  by 
Louis  et  al.  (1988)  and  slightly  different  models  and  parameters  for  the 
satellite  height  observational  errors  were  used  in  the  OSEs  reported  here. 
These  observational  errors  are  described  in  detail  in  Section  2.3  of  Louis  et 
al.  (1988). 
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In  the  SSM/T  OSSEs  we  do  not  directly  use  the  statistical  properties 
described  in  Section  5  which  were  used  to  simulate  the  data.  Instead  we  saved 
the  actual  errors  used  in  simulating  the  data  and  used  them  to  develop  global 
models  of  the  observing  errors  for  use  in  the  01.  Our  first  motivation  for 
thi.s  is  that  the  results  of  Section  5  do  not  allow  us  to  estimate  horizontal 
correlations.  Secondly,  the  models  of  the  observing  errors  used  in  our  sensor 
simulation  study  are  considerably  more  complex  than  those  used  by  the  ASAP  01 
and  we  did  not  wish  to  make  major  modifications  to  the  analysis  procedure.  In 
studying  the  global  simulated  errors  we  grouped  all  errors  generated  for  the 
10  and  12  GMT  ir.tervals  together.  Basically  we  found  that; 

1)  The  vertical  height  error  correlations  were  fairly  well  fit  by  a  simple 
func  t ion . 

V  -  1/(1  +  k(A(ln  p)!^ 

where  we  found  k  to  be  .3744  by  a  least  squares  fit.  The  observed 
correlations  and  the  fitted  function  are  shown  in  Fig.  6.1.  The 
corresponding  RH  correlation  were  all  small  and  deemed  not  significantly 
different  from  zero. 

2)  The  horizontal  height  error  correlation  are  close  to  zero.  The  RH 
horizontal  error  correlations  are  also  very  small  except  at  300  and 
400  mb  where  they  remain  above  .2  out  to  3000  km  (Fig.  6.2).  In  both 
cases  we  modeled  the  horizontal  error  correlations  by 

n  -=  exp(  -  (d/dg)^) 

where  we  choose  dg  as  83  km  for  height  and  198  km  for  RH. 

3)  No  significant  biases  were  found  and  the  rms  errors  were  of  the  expected 
size  (see  Table  6.1). 

The  .statistical  models  used  in  the  01  make  use  of  these  findings. 

The  ASAP  01  was  adjusted  to  assimilat  the  SSM/T  data.  The  SSM 
retrieval-s  of  RH  on  mandatory  pressure  levels  are  used  directly,  while  the 
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T  retrievals  are  first  converted  to  heights.  In  principle,  this  transforma¬ 
tion  is  not  necessary,  the  analysis  scheme  could  make  use  of  the  the  fact  that 
sigma  level  heights,  the  analyzed  variables,  are  correlated  with  the  observed 
mandatory  level  or  layer  temperatures.  For  SSM/T  we  retrieve  temperature  on 
mandatory  levels  directly,  instead  of  layer  temperatures.  The  level  tempera¬ 
tures  are  interpolated  linearly  in  In  p  to  the  sigma  levels  and  then  inte¬ 
grated  hydrostatically  to  obtain  the  sigma  level  heights.  The  predicted  or 
first  guess  height  at  the  first  sigma  level  within  the  retrieved  temperature 
profile  is  used  to  begin  the  integration.  The  Integration  also  assumes 
temperature  varies  linearly  in  In  p.  The  procedure  used  is  identical  to  the 
procedure  for  anchoring  satellite  thicknesses  described  in  Hoffman  et  al . 
(1988)  from  the  point  at  which  mandatory  level  temperatures  are  available. 

Note  that  the  retrieved  profile  is  not  extrapolated  vertically  at  all. 

The  data  selection  algorithm  was  altered  to  recognize  the  SSM/T  data. 

The  basic  data  selection  algorithm  follows  Bergman  (1979)  as  described  by 
Norquist  (1988).  The  first  stage  of  data  selection  is  done  in  terms  of 
profiles:  Data  items  in  the  up  to  8  closest  profiles  are  candidates  to  be 
chosen  for  use  in  the  actual  analysis  in  the  second  stage.  Here  closeness  is 
measured  by  the  magnitude  of  height  height  forecast  error  correlation  times 
the  number  of  non  missing  data  items  in  the  profile  times  a  measure  of  data 
quality.  (A  Z  observation  is  considered  one  data  item  and  a  (u,v)  pair  is 
considered  one  data  item.)  The  data  quality  is  taken  to  be  unity  for  all  data 
except  it  is  0.42  for  TOVS  and  0.92  for  SSM/T  height  data.  These  data  quality 
values  are  rough  estimates  of  the  ratio  of  RAOB  to  TOVS  and  RAOB  to  SSM/T 
height  OESDs .  Thus  RAOB  and  other  type  1  data  are  preferred  over  all 
others.  A  SSM/T  profile  will  be  preferred  to  a  TOVS  profile.  The  second 
stage  of  data  selection,  which  was  not  altered,  selects  for  each  analysis 
point  up  to  10  data  items  from  the  selected  profiles  which  individually  would 
give  the  largest  reduction  in  estimated  analysis  error.  In  this  second  stage 
Che  actual  estimated  OESDs  are  used. 

6.1.2  Forecast  and  initialization 

The  AFGL  normal  mode  initialization  (NMI)  is  based  <-«n  the  NMC  NMI 
(Ballish,  1980).  The  AFGL  global  spectral  model  is  based  on  the  NMC  GSM 
designed  by  Sela  (1980).  For  the  version  used  here,  the  physics  routines  are 
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taken  almost  intact  from  NMC  (circa  1983).  The  hydrodynamics,  i.e.,  the 
adiabatic,  inviscid  dynamics  including  vertical  and  horizontal  advectlon,  time 
stepping,  and  transformations  between  spectral  and  physical  space,  were 
completely  redesigned,  as  documented  by  Brenner  et  al.  (1982,  1984). 

Tiiere  are  a  number  of  parameters  in  the  forecast  and  initialization  codes 
that  can  be  adjusted.  Briefly,  the  spectral  resolution  of  the  forecast  model 
is  defined  by  a  rhomboidal  truncation  at  wave  number  30.  The  Gaussian  grid  of 
the  forecast  model  contains  76  x  96  latitude  longitude  points.  There  are  12 
layers,  the  first  (top)  5  of  which  have  no  moisture.  The  sigma  interfaces  are 
at  O.CO,  0.05,  0.10,  0.15,  0.20,  0.25,  0.30,  0.375,  0.50,  0.65,  0.80,  0.925, 
and  1.00.  The  time  scheme  used  is  centered  semi  -  implicit  with  a  time  step  of 
17.25  minutes.  Horizontal  fourth  order  diffusion  (-kV^)  is  applied  to  all 
modes  of  divergence  and  to  modes  in  the  upper  half  of  the  rhomboid  for 
vorticity,  temperature  and  specific  humidity.  The  diffusion  coefficient  used 
is  K  -  6.10^^  ra^  s’^.  In  the  NMI ,  two  Machenauer  iterations  are  applied  to 
modes  for  the  four  largest  equivalent  depths  which  have  periods  less  than  or 
equal  to  48  h. 

6.2  Verification  Procedures 

Our  verification  procedures  include  subjective  comparisons  of  analyses 
and  forecasts,  quantitative  comparisons  of  rms  errors  of  analyzed  and  forecast 
fields  and  a  calibration  of  the  quantitative  measures  making  use  of  the  OSE 
results.  These  procedures  will  be  described  in  more  detail  in  the  next 
section.  Here  we  describe  the  diagnostic  fields  which  we  compare. 

One  of  the  performance  measures  we  use  in  measuring  impact  is  the  rms  and 
bias  error  of  a  diagnosed  total  cloud  fraction.  Although  the  "real"  simulated 
clouds  generated  during  the  nature  forecasts  were  not  saved,  the  large  scale 
(i.e.,  nonconvective)  cloud  fraction  is  readily  diagnosed  from  the  RH  field. 
This  sort  of  measure  puts  a  premium  on  proper  foreca.sting  at  the  higher  range 
of  RH,  since  errors  at  lower  values  have  no  effect  on  large  scale  cloud 
amounts . 

To  convert  RH  to  level  cloud  coverage  fraction,  f ^ ,  we  invert  the  Tibaldi 
scheme  described  by  Norquist  (1988,  Appendix  A).  Level  clouds  are  obtained  at 
the  first  six  mandatory  layers,  i.e.  1000,  850,  700,  500,  400,  and  300  mb. 
There  are  combined  pairwise  assuming  maximal  overlap  into  low,  middle  and  high 
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layer  cloud  fractions,  and  fjj.  The  three  layer  cloud  fractions  are 

then  combined  into  a  total  cloud  fraction  fj  assuming  random  overlap, 

(1  -  f-^)  -  (1  -  f^Xl  -  £„)<!  -  fH) 

6.3  Spinup  experiment 

To  begin  our  experiments  we  first  performed  two  96  h  forecasts  starting 
from  "perfect"  initial  conditions  on  00  GMT  11  and  21  November  1979.  These 
initial  conditions  are  taken  directly  from  the  nature  run,  interpolated  to  the 
model  sigma  structure,  analyzed  into  spherical  harmonic  coefficients  and 
initialized  with  the  adiabatic  NMI .  We  found  the  growth  of  forecast  error 
(i.e.,  forecast  -  nature)  to  be  rather  slow  in  these  forecasts.  In  fact  the 
forecast  ending  on  the  00  GMT  15  November  was  not  much  worse  than  a  typical 
analysis.  Accordingly,  this  state  was  used  as  the  initial  "analysis"  for  a 
three  day  SPINUP  assimilation  experiment,  ending  00  GMT  18  November.  SPINUP 
makes  use  of  the  standard  STATSAT  configuration. 

6.4  Schedule  of  impact  experiments 

The  end  of  the  SPINUP  assimilation  is  used  as  the  starting  analysis  for 
all  our  OSSEs  which  therefore  run  from  00  GMT  18  November  through  00  GMT 
25  November.  (The  first  analysis  of  each  OSSE  is  at  06  GMT  18  November.)  For 
each  OSSE  96  h  forecasts  are  made  from  00  GMT  21,  23  and  25  November. 

The  OSSEs  described  here  are  STATSAT,  NCSAT,  STATSAT  +  SSM/T  and  NOSAT  + 
SSM/T.  STATSAT  includes  all  the  Level  II  data  which  were  simulated  by  NMC  as 
described  in  Section  3.2  except  that  surface  observations  are  not  used  and 
satellite  temperature  soundings  over  land  are  not  used.  In  NOSAT  the 
satellite  temperature  soundings  and  CDW  observations  are  excluded.  In  the 
other  experiments  the  SSM/T- 1  and  T-2  data  are  added. 

The  OSEs  STATSAT  and  NOSAT  are  analogous  to  the  OSSEs,  at  least  as  far  as 
data  usage  is  concerned.  The  OSEs  however  were  run  for  one  week  each  during 
February  and  June  1979.  These  experiments  are  described  in  detail  by  Louis  et 
al .  (1988)  and  are  used  here  primarily  to  calibrate  the  OSSE  results.  One 
notable  differenc  .  between  the  OSSE  and  OSE  experiments  is  that  the  OESD  for 
satellite  heights  for  STATSAT  and  NOSAT  in  the  February  OSE  were  substantially 
larger  than  in  the  other  experiments.  (These  values  are  given  by  Louis  et 
al.,  1988.) 
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Table  6.1 


SSN/T  OESDs  used  in  data  assimilation  experiments 


p  ( mb ) 

1000 
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700 

500 

400 

300 

250 

200 

150 

100 

70 

50 

Z(m) 

7 . 98 

7.98 

14.01 

21.71 

25.08 

26.57 

27.01 

28.38 

32.15 

38.37 

39.40 

41.85 

RH(%}  7.91 


18.35  174.45  15.12  15.17  19.35 
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Figure  6.1  Vertical  height  error  correlations  of  simulated  SSM  retrievals 
based  on  collocation  with  radiosondes.  Solid  curve  Is  least 
squares  fit  to  data  used  In  01. 
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Figure  6.2  Horizontal  height  error  correlations  of  simulated  SSM  retrievals 
based  on  collocation  with  radiosondes.  Shown  are  observed  errors 
at  300  mb,  500  mb  and  function  used  in  01. 
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7 .  OSSE  Results 


7.1  Subjective  synoptic  evaluation 

The  spinup  forecast,  starting  from  perfect  initial  conditions  at  00  GMT 
11  November,  and  the  spinup  assimilation,  which  was  performed  from  the  end  of 
the  spinup  forecast  at  00  GMT  15  November  to  00  GMT, 18  November  were  compared 
with  the  corresponding  nature  data.  The  spinup  forecast  had  very  similar 
500  mb  height  patterns,  but  our  forecast  Is  consistently  warmer  than  the  ECMWF 
forecast.  The  possible  reason  for  this  warm  bias  is  the  lack  of  a  radiation 
parameterization  in  our  model;  a  typical  radiative  cooling  of  the  atmosphere 
of  1  K/day  would  correspond  to  a  500  mb  height  difference  of  roughly  20  m/day, 
which  is  consistent  with  the  approximate  height  difference  of  80  m  at  the  end 
of  the  4-day  forecast.  The  analyses  during  the  spinup  assimilation,  which 
correspond  to  a  STATSAT  configuration,  are  also  quite  similar  to  the  nature 
data,  except  that  they  are  considerably  noisier,  possibly  due  to  the  warm  bias 
of  the  first  guess,  which  is  only  corrected  at  data  locations. 

Results  for  the  Northern  Hemisphere  from  the  OSSE  assimilation  period 
(00  GMT  18  November  to  00  GMT  25  November)  are  shovm  in  Fig.  7.1  -  7.4  for  00 
GMT  23  November.  The  500  mb  height  pattern  of  the  nature  data  (Fig.’  7,1a) 
shows  a  distinct  wavenumber  four  pattern,  which  is  present  throughout  the 
entire  November  time  period.  Several  smaller  scale,  mobile  troughs  are 
superimposed  on  the  long-wave  structure.  The  analyses  from  all  three  OSSEs 
shown  in  Fig.  7.2  -  7.4  are  visibly  noisier  than  nature,  both  at  1000  and 
500  mb.  The  error  patterns  of  all  three  OSSEs  look  similar:  the  long-wave 
low  over  Eastern  Canada  is  too  low  in  the  OSSEs,  and  there  is  large  area  of 
negative  height  errors  at  70°-80°N,  centered  at  90’E.  These  two  features  are 
present  throughout  the  entire  assimilation  period.  The  error  over  Canada  is 
worst  in  STATSAT,  whereas  the  other  feature  is  worse  in  the  SSM  OSSEs.  There 
are  two  regions  where  the  SSM+TOVS  analysis  is  noticeably  worse  than  either 
SSMSAT  or  STATSAT:  one  is  the  long-wave  trough  at  150*W,  which  is  much  too 
deep  in  SSM+TOVS,  less  so  in  SSMSAT,  and  approximately  correct  in  STATSAT,  the 
other  a  small  short  wave  over  Newfoundland,  which  is  not  captured  well  in  any 
of  the  OSSEs,  with  the  largest  errors  in  SSM+TOVS,  smaller  errors  in  STATSAT, 
and  the  smallest  errors  in  SSMSAT.  At  day  7  of  the  assimilation  (00  GMT  25 
November,  Fig.  7.5  -  7.8)  the  500  mb  height  errors  over  the  Pacific  have 
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grown,  and  tlioy  now  appear  in  all  three  OSSEs.  Over  the  Eastern  Atlantic  is 
ar.  ar.wi  of  large  height  errors  in  SSM+TOVS  (Fig.  7.8),  which  is  virtually 
ab.sent  in  either  SSMSAT  (Fig.  7.7)  or  STATSAT  (Fig.  7.6).  The  fact  that 
adding  TOVS  data  to  Che  SSMSAT  data  lead  to  such  a  drastic  degradation  of  the 
ar.ilyst-s  in  those  regions  is  troubling,  since  the  01  is  designed  to  optimally 
combine  all  available  data,  given  the  correct  statistics  of  the  true  field.s 
and  errors. 

Some  possible  reasons  for  this  analysis  degradation  are  related  to  the 
fact  that  the  satellite  thicknesses  are  anehorod  to  the  first  gues.s  surface 
pre.ss'are.  Satellite  data  are  thus  not  able  to  correct  barotropic  first  guess 
height  errors,  i.e.  errors  present  at  both  the  1000  and  500  mb  levels.  The 
error  stati.stics  used  in  the  01,  however,  do  not  take  this  anchoring  error 
irito  account,  thus  resulting  in  inappropriately  large  weights  being  given  tn 
.^.itollite  "height"  ooservations  ,  at  the  expense,  of  other,  non  -  sate  1 1  i  te 
data.  Even  in  the  absence  of  nonsatellite  data,  the  addition  of  satellite 
data  may  worsen  the  500  nab  height  analysis,  if  the  first  guess  errors  in 
surface  pressure  and  500-1000  mb  thickness  are  compensating  each  other.  All 
these  potential  problems  would  be  exacerbated  by  adding  TOVS  to  the  SSM 
data.  Another  possibility  is  that  in  the  presence  of  significant  biases  in 
the  satellite  data,  which  we  have  ignored  in  our  OSSEs,  mixing  TOVS  and  SSM 
data  in  overlapping,  data  swaths  may  introduce  noise  in  the  analysis.  Finally, 
because  of  central  memory  limitations,  some  aircraft  and  cloud  track  winds 
used  in  the  SSMSAT  analyses  could  not  be  used  in  SSM+TOVS.  Any  of  these 
reasons  may  result  in  initially  small  analysis  differences,  which  could  be 
amplified  during  the  6 -hour  forecast  between  analysis  time  periods,  and 
through  differences  in  data  selection  due  to  quality  control  procedures  that 
depend  on  the  first  guess.  In  an  effort  to  identify  these  scenarios,  analyses 
and  first  guess  fields  for  the  OSSEs  were  studied  along  with  the  corresponding 
nature  data  and  the  simulated  observation  residuals. 

The  analysis  error  over  the  Atlantic  and  the  short  wave  associated  with 
it,  which  at  00  GMT  25  November  is  lov.ated  over  Ireland  in  the  nature  data, 
but  is  displaced  about  10°  West  in  the  SSM+TOVS  analysis,  can  be  traced  back 
to  Newfoundland  at  00  GMT  23  November,  at  which  point  the  analysis  errors  for 
the  SSM+TOVS  analyses  are  already  substantially  larger  than  for  the  SSMSAT 
analyses.  The  two  analyses  diverge  over  that  area  in  the  preceding  26 
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hours:  the  500  mb  height  errors  at  00  GMT  22  November  (Fig.  7.9)  are  much 

more  similar  in  SSMSAT  and  SSM+TOVS  over  Newfoundland.  In  the  time  period 
between  00  GMT  22  November  to  00  GMT  23  November,  SSM  and  TOVS  data  are  both 
present  over  the  area  only  at  06  GMT  and  18  GMT,  i.e.  at  times  when  there  are 
no  RAOBs  over  the  area;  the  impact  of  the  TOVS  data  at  those  times  should  be 
small,  given  the  higher  data  density  and  smaller  OESDs  of  the  SSM  data.  At 
the  other  times,  only  TOVS  satellite  data  are  present  over  the  oceans,  but 
there  were  no  obvious  instances  where  the  additional  thickness  information 
might  have  lead  to  larger  500  mb  height  errors.  It  is  thus  not  entirely  clear 
what  caused  the  two  analyses  to  diverge  over  that  area. 

Whatever  the  reasons  for  initial,  small  analy^j-S  differences  are,  the 
persistence  and  amplification  of  these  errors  requires  some  additional 
explanation.  As  is  shown  in  Fig.  7.8b,d  and  7.10a,b,  the  analysis  errors  and 
the  first  guess  errors  at  00  GMT  23  November  have  a  large  barotropic 
component,  and,  as  a  consequence,  the  analysis  is  ineffective  over  the  ocean 
areas,  where  there  are  no  radiosonde  observations  of  geopotential.  In 
agreement  with  this  scenario,  the  satellite  observation  residuals  are  small 
compared  to  the  actual  analysis  errors.  Where  radiosondes  are  present,  the 
analysis  errors  are  smaller  than  the  first  guess  errors  right  over  those 
locations,  except  one  radiosonde  at  12  GMT  22  November  (at  Sept  Isles,  Quebec, 
at  50.22‘’N,  66.25'’W)  and  another  at  00  GMT  23  November  (St.  Johns, 
Newfoundland,  at  A7.62°N,  52.75°W),  which  were  rejected  by  the  gross  error 
check  in  SSM+TOVS,  but  not  in  SSMSAT.  At  those  locations,  the  first  guess 
field  in  SSM+TOVS  had  diverged  too  far  from  the  truth  for  the  RAOBs  to  be 
used.  The  rapid  error  growth  during  this  short  period  is  taking  place  in  a 
strongly  baroclinic  zone,  in  which  a  strong  sensitivity  to  initial  conditions 
of  even  a  ur  forecast  can  be  expected.  After  00  GMT  23  November,  the 
short  wave  is  entirely  over  the  ocean,  thus  making  it  impossible  for  the 
analysis  to  correct  the  by  now  well  established  height  errors. 

The  other  area  of  large  analysis  errors  apparent  in  Fig.  7.8b  is  over  the 
Pacific,  near  the  international  dateline.  This  error,  which  is  associated 
with  a  longwave  trough  in  the  nature  data,  can  also  be  traced  back  several 
davs  (it  is  already  visible  at  00  GMT  21  November),  but  it  is  an  essentially 
stationary  feature.  The  longwave  trough  Is  consistently  analyzed  too  deeply 
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(in  bor'i  SSMSAT  and  SSMtTOVS  analyses),  and  at  several  times  (including  00  GMT 
23  Noveraber)  a  smaller  scale  ridge  downstream  of  the  trough  is  too  strong. 

The  troui.h  is  always  analyzed  more  deeply  in  SSM+TOVS  than  SSMSAT.  As  was  the 
case  over  the  Atlantic,  the  analysis  errors  have  a  strong  barotropic  com¬ 
ponent  Sinte  the  longwave  trough  is  entirely  over  the  ocean,  only  satellite 
height  data  are  available,  and  again  the  height  residuals  were  found  to  be 
much  smaller  than  tlte  actual  first  guess  height  errors  because  of  the  anchor¬ 
ing  of  the  satellite  th.i  cknesses .  A  comparison  of  the  first  guess  a-ij  analy¬ 
sis  errors  at  00  GMT  23  .November  show  the  negligible  impact  of  the  analysis, 
ilIowiTAg  the  persistence  of  the  first  guess  errors.  The  reason  for  the  ini¬ 
tial  d..verger.i'.e  of  the  t’S.MSA'I  and  SSM+TOVS  analyses  is  again  not  entirely 

i  t*  I  . 

.Ai'i  exarpl.:-  of  the  Southern  Hemisphere  height  fields  is  given  in 
lir, .  7.11  -  7,14  for  00  GMT  23  November.  Again,  the  OSSF.  analyses  arc  all 
ir.ucli  noisier  than  nature.  In  addition,  the  500  mb  lioights  are  consistently 
too  high  over  Antarctica.  The  nature  data  show  a  number  of  small  scale, 
mobile  troughs  at  300  and  1000  mb,  a  situation  typical  for  the  entire 
assimilation  period.  Analysis  errors  are  due  to  both  amplitude  and  phase 
errors  of  the.s  ■;  features,  and  are  present  at  both  500  and  1000  mb.  .Overall, 
botii  the  SSM+TOVS  an.i  SSMSAT  analyses  seem  the  best,  and  the  STATSAT  analyses 
worst  over  the  Soutf:crn  Hemisphere. 

As  discussed  in  previous  sections  the  OSSE  included  simulated  retrievals 
of  relative  humidity  from  SSM/T-2  brightness  temperatures.  Therefore,  it  is 
also  of  interest  to  assess  the  impact  of  such  data  upon  the  analyses  and 
forecasts  of  moisture  during  the  assimilation  period.  Comparisons  between  the 
analyses  of  different  experiments  is  not  very  informative  since  instantaneous 
fields  of  relative  humidity  are  rather  noisy.  Instead,  we  display  only 
difference  plots,  averaged  fields,  and  differences  of  averaged  fields  for 
illustrative  purposes. 

In  Figs.  7.15  and  7.16  we  show  difference  plots  of  the  analyzed  moisture 
field  at  850  mb  with  the  nature  run  field  for  2  different  times  during  the 
assimilation  period  (23  November  and  25  November,  both  at  00  GMT).  In  each 
figure  fields  for  the  STATSAT,  SSMSAT,  and  SSM+TOVS  experiments  are  shown. 

In  this  set  of  figures  differences  in  relative  humidity  are  plotted  with 
a  contour  interval  of  25%  with  negative  differences  shovm  as  dashed  line.s  am) 
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positive  differences  shovm  as  solid.  The  times  shovm  represent  analysis 
differences  at  days  5  and  7  of  the  assimilation  period.  It  is  quite  clear 
that  the  moisture  analyses  created  with  the  SSM  retrievals  are  much  improved 
over  the  STATSAT  analyses.  In  addition,  the  results  for  SSMSAT  and  SSM+TOVS 
are  very  similar  with  only  minor  differences  seen  between  the  two  experi¬ 
ments.  It  seems  that  the  presence  or  absence  of  TOVS  retrievals  has  only  a 
small  impact  upon  the  humidity  analysis  and  that  nearly  all  the  improvement 
stems  from  inclusion  of  the  additional  SSM  data.  Improvement  is  most  dramatic 
over  southern  hemisphere  ocean  areas  where  conventional  moisture  data  is 
practically  nonexistent.  Smaller,  but  equally  clear  improvements  are  seen  in 
the  northern  hemisphere  as  well.  For  example,  note  the  improved  analysis  over 
the  central  North  Atlantic  Ocean  on  25  November.  Also  noteworthy  is  the  fact 
chat  the  analyses  are  improved  over  land  areas  even  though  SSM  retrievals  were 
confined  to  Che  oceans.  This  is  quite  evident  over  equatorial  Africa  on 
25  November,  and  over  South  America  on  23  November.  Thus  the  use  of  the 
additional  data  source  over  one  region  has  the  potential  to  improve  subsequent 
analyses  in  more  distant  areas  downstream  from  the  data  location. 

Another  useful  means  of  viewing  the  moisture  fields  associated  with  the 
OSSEs  is  by  converting  relative  humidity  fields  to  cloud  cover  fields.  The 
resulting  cloud  cover  fields  offer  some  insight  into  the  accuracy  of  the 
analysis/forecast  relative  humidity  fields  and  the  algorithm  which  is  used  for 
the  humidity  to  cloud  conversion. 

We  have  used  moisture  fields  from  both  the  STATSAT  and  SSMSAT  experiments 
to  infer  large  scale  cloudiness.  For  this  work  the  inverse  Tibaldi  scheme  was 
applied  (see  Section  6,2).  We  found  the  complete  relative  humidity  fields 
(i.e.  with  all  spectral  coefficients  retained)  to  be  highly  noisy  making 
interpretation  difficult.  Instead,  we  filtered  the  fields  by  using  a  T15 
spectral  truncation  which  allows  us  to  focus  only  on  the  large  scale  features. 

Cloud  cover  results  for  STATSAT  and  SSMSAT  at  00  GMT  25  November  are 
shown  in  Figures  7.17  and  7.18,  respectively.  We  show  only  the  inferred  low 
cloud  cover  with  a  contour  interval  of  25  percent.  For  a  surface  pressure  of 
1000  mb  this  would  include  clouds  in  the  layer  from  950  mb  to  733  mb.  In 
Fig.  7.17  is  the  cloud  cover  inferred  from  the  nature  run.  Remaining  panels 
show  the  inferred  cloudiness  of  the  analyses  valid  at  the  same  time  and  the 
corresponding  analysis  -  nature  difference  fields. 
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Main  features  in  the  nature  run  cloudiness  include  tha  low  amounts  ot 
cloudiness  over  the  northern  hemisphere  subtropics  associated  with  the 
desccTidlng  branch  of  the  Hadley  circulation.  Note  in  particular  regions  with 
less  than  25%  cloud  cover  over  the  southwest  U.S.,  northern  Africa,  and 
southern  Asia.  High  cloud  amounts  aiv.  seen  in  the  equatorial  regions  of 
Africa  as  well  as  South  America  which  It  might  be  noted,  are  close  to  the 
general  location  of  the  ITCZ  as  it  advances  southward  during  November. 

As  seen  in  the  relative  humidity  fields,  the  SSMSAT  results  are  more 
realistic  (relative  to  the  nature  run)  than  those  of  STATSAT.  In  SSMSAT  the 
analy.si.s  of  cloudiness  over  central  Africa  and  South  America  is  in  reasonable 
agreement  with  the  nature  run  while  the  STATSAT  field  completely  misses  the 
high  cloud  amounts  in  these  same  areas. 

The  difference  plots  also  show  the  SSMSAT  analysis  to  be  better  than 
STATSAT.  Again  note  South  America  and  equatorial  Africa  where  the  cloud 
amount  errors  in  STATSAT  exceed  -75%  over  large  areas.  It  is  also  interesting 
to  note  that  nearly  all  errors  in  STATSAT  and  SSMSAT  over  the  continents  are 
negative.  Tnis  is  most  likely  a  reflection  of  the  fact  that  the  physics 
package  in  the  GSM  contains  no  evaporation  over  land  surfaces;  consequently 
the  lowest  levels  in  the  model  tend  to  be  too  dry. 

Finally,  although  we  found  that  the  inferred  cloud  coverage  agrees  with 
the  large  scale  relative  humidity  fields  in  a  qualitative  sense,  we  also  £o  .ind 
that  cloud  amounts  were  generally  much  too  high.  Global  statistics  for  the 
nature  run,  assimilations,  and  forecasts  showed  an  average  total  cloud  cover 
of  roughly  80  percent  while  climatological  values  are  known  to  be  closer  to 
50  percent.  To  the  extent  that  that  the  resulting  cloud  fields  mirror 
relative  humidity  patterns  this  implies  that  humidity  in  both  the  GSM  and 
nature  run  is  too  high.  However  it  is  also  possible  that  the  humidity  to 
cloud  conversion  method  is  biased.  Any  quantitative  use  of  the  cloudiness 
inferred  by  this  means,  for  example  simulation  of  satellite  imager  data,  will 
require  some  tuning  of  the  scheme  to  reduce  the  observed  bias. 
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7.2  Objective  statistical  evaluation 

As  an  objective  measure  of  the  quality  of  the  analyses  and  forecasts, 
root  mean  square  (rms)  and  mean  (bias)  statistics  of  the  differences  from 
nature  were  calculated.  Fig.  7.19  shows  the  global  500  mb  rms  height  error  of 
the  NOSAT,  STATSAT,  SSM,  and  SSM+TOVS  analyses  and  forecasts.  The  NOSAT 
analysis  errors  increase  from  Che  35  m  typical  for  the  STATSAT  analysis  to 
50  m  by  day  4  of  the  assimilation,  whereas  the  STATSAT  analysis  errors 
decrease  by  only  1-2  m  over  the  assimilation  period.  The  forecast  error 
growth  is  more  rapid  in  STATSAT,  but  errors  remain  smaller  than  those  of  the 
NOSAT  forecast  for  the  entire  4-day  forecasts.  The  SSM  analysis  errors  are 
consistently  smaller  than  those  of  STATSAT,  by  up  to  2  to  3m.  The  day  3 
forecast  errors  are  smaller  than  those  of  STATSAT,  but  the  day  7  forecast  is 
worse.  The  SSM+TOVS  analysis  have  larger  errors  than  those  of  SSM,  and  the 
forecasts  are  either  the  same  (day  7)  or  worse  (day  3) .  Compared  to  STATSAT, 
SSM+TOVS  analyses  and  the  day  3  forecast  is  only  slightly  worse  (by  less  than 
1  m) ,  whereas  the  day  7  forecast  is  noticeably  worse.  The  objective  results 
thus  confirm  our  impression  from  the  subjective  evaluation  that  the  SSMSAT 
analyses  (and  forecasts)  are  closest  to  nature.  The  reasons  for  the 
degradation  of  the  analyses  in  SSM+TOVS  discussed  in  the  previous  section  in  a 
regional  context  thus  apply  globally,  as  well. 

Results  at  other  levels  largely  mirror  those  at  500  mb.  The  1000  mb 
height  statistics  (not  shown)  show  a  much  smaller  impact  of  the  satellite  data 
(NOSAT  analysis  errors  differ  by  no  more  than  3  ra  from  STATSAT),  but 
qualitatively  the  same  results  apply.  In  particular,  the  SSMSAT  analyses  are 
closer  to  nature  than  either  STATSAT  or  SSM+TOVS,  the  day  3  forecast  is  most 
skillful  in  SSMSAT,  whereas  the  day  7  forecast  is  roughly  the  same. 

The  global  error  statistics  of  humidity  analyses  and  forecasts  for 
STATSAT,  SSMSAT,  and  SSM+TOVS  are  shown  in  Fig.  7.20.  The  rms  error  curves 
for  all  three  experiments  were  obtained  by  differencing  the  850  mb  analyses 
and  forecasts  with  the  corresponding  nature  run  forecasts.  In  each  plot  the 
solid  curves  denote  the  assimilation  errors  and  the  broken  curves  the  errors 
of  the  forecasts  which  were  run  off  the  day  3,  5,  and  7  analyses. 

The  error  stati.stics  for  the  STATSAT  and  NOSAT  experiments  may  be 
considered  baselines  against  which  to  judge  the  other  OSSEs.  In  STATSAT 
simulated  TOVS  temperature  retrievals  were  included  with  the  rest  of  the 
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conventional  radiosonde  network  data  used  in  the  NOSAT  experiment.  Since  the 
STaTSAT  experiment  used  the  same  moisture  data  we  would  expect  the  relative 
humidity  error  statistics  to  be  quite  similar.  In  fact,  comparisons  of 
STATSAT  with  NOSAT  (not  shown)  show  the  analysis  errors  to  be  virtually 
identical,  with  differences  in  the  rms  errors  of  less  than  1%.  The  actual 
values  of  the  rms  errors  throughout  the  7  day  assimilation  period  are 
approximately  27%.  Errors  in  the  forecasts  approach  30%  by  4  days. 

In  the  SSMSAT  experiment  the  simulated  SSM  temperature  and  moisture 

retrievals  were  inserted  at  the  same  time  that  the  TOVS  data  was  excluded. 

The  results  showt'  clearly  indicate  an  improvement  in  the  accuracy  of  the 
moisture  analyses  and  forecasts  when  such  data  is  used.  Analysis  errors  are 
roughly  5%  smaller  throughout  the  7  day  period  with  rms  errors  values  of  about 

22%  Forecast  errors  are  also  initially  less  than  those  of  STATSAT  by  5%  but 

increase  more  quickly  so  that  after  4  days  they  differ  by  only  2%. 

In  SSM+TOVS  both  types  of  satellite  data  were  included.  As  expected 
based  on  the  subjective  evaluations  the  Impact  of  adding  the  TOVS  data  to  the 
SSM  retrievals  is  at  beat  minimal.  Both  analyses  and  forecasts  for  the 
SS.M+TOVS  experiment  have  slightly  larger  rms  errors,  but  in  all  instances  the 
differences  are  less  than  or  equal  to  2%.  Note  the  tendency  for  the'  forecast 
error  curves  to  flatten  out  late  in  the  forecast  period  at  about  30%.  This 
most  likely  represents  the  upper  limit  on  possible  analysis/forecast  errors 
based  on  the  natural  variability  of  the  relative  humidity  fields. 

In  sum,  both  the  subjective  and  objective  evaluations  of  the  moisture 
data  indicate  a  distinct  positive  impact  upon  relative  humidity  analyses  when 
simulated  SSM-derived  moisture  retrievals  are  added  to  conventional 
observations  during  the  global  assimilation  cycle.  Only  minor  impact 
(relative  to  SSMSAT)  is  seen  when  the  SSM  data  is  combined  with  TOVS 
temperature  retrievals. 

7.3  Evaluation  of  zonal  cross  sections 

We  examined  zonal  cross  sections  of  u  and  v  wind  components,  temperature 
and  relative  humidity  at  individual  synoptic  times  and  averaged  over  the  last 
five  0000  GMT  analyses  of  the  experiments.  Zonal  averaging  is  denoted  hero  by 
square  bracKets  ([  ])  and  time  averaging  by  an  overbar  (  ).  Our  purpose  hero 
is  to  determine  how  well  the  assimilation  system  is  capturing  the  mean  meri - 
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dional  circulation,  pole  to  equator  temperature  and  humidity  structure  and 
zonal  jets  and  to  describe  the  Impact  of  the  different  observing  systems  on 
these  features.  These  features  are  important  climate  diagnostics.  Errors  in 
zonally  averaged  quantities  are  also  important  to  NWP  because  errors  in  the 
climate  make  a  persistent  contribution  to  errors  In  the  analyses  and  these 
errors  may  be  useful  in  diagnosing  faults  in  the  physical  parameter izations 
used  in  the  model.  In  the  following,  the  impact  of  the  SSM  OSSEs  will  be 
compared  with  that  of  the  WINDSAT  OSSE.  The  WINDSAT  OSSE  is  identical  to 
STATSAT,  except  that  simulated  wind  observations  from  a  proposed  satellite- 
based  lidar  wind  profiler  were  made  available  to  the  assimilation.  The 
WINDSAT  OSSE  is  described  in  more  detail  in  a  separate  technical  report 
(Hoffman  et  al.,  1989).  We  concentrate  here  on  the  time  averaged  fields  for 
t.I._  nature  run  and  SiAtbAi  and  on  tne  impact  of  SSM  on  the  wind,  temperature 
and  humidity  fields.  These  fields  and  differences  from  the  nature  run  are 
displayed  in  Figs.  7.21  through  7.24.  For  comparison  the  some  results  for  the 
WINDSAT  experiment  are  also  displayed. 

Considering  first  the  zonal  time  averaged  zonal  wind  component 
(Fig.  7.21),  we  see  in  all  cases  the  midlatitude  jets  peak  near  200  mb.  The 
N.H.  jet  is  somewhat  narrower  and  stronger.  The  jet  maximum  is  neat  30  m/s  in 
all  cases.  Easterly  winds  extend  through  the  depth  of  the  atmosphere  in  the 
tropics.  The  trade  winds  (surface  easterlies)  are  a  bit  stronger  in  the 
N.H.  The  nature  [u]  is  quite  similar  to  observed  fields.  Compared  to  the 
GFDL  monthly  mean  [u]  for  November  1979  (Lau,  1984),  the  jets  in  the  nature 
run  are  positioned  somewhat  poleward  and  have  slightly  different  magnitudes. 
STATSAT,  SSMSAT  and  WINDSAT  [u]  agree  well  with  the  nature  run.  In  all  three 
cases  the  errors  are  order  1  m/s  with  WINDSAT  having  somewhat  smaller  errors. 

In  Fig.  7.22,  which  depicts  [v]  both  Hadley  and  Ferrel  cells  are  readily 
apparent.  On  the  other  hand  the  surface  southerlies  in  the  S.H.  Ferrel  cell 
are  stronger  than  in  the  N.H.  The  maximum  convergence  at  the  surface  where 
the  two  returning  branches  of  the  Hadley  cell  meet  is  at  10°N.  The  cor¬ 
responding  upper  level  divergence  is  at  250  mb.  The  N.H.  Hadley  cell  appears 
to  be  somewhat  stronger  than  the  S.H.  Hadley  cell.  Compared  to  the  GFDL 
monthly  means,  the  surface  winds  agree  fairly  well,  but  the  poleward  branches 
of  the  Hadley  cell  are  twice  as  strong  in  the  GFDL  analyses.  Comparing 
STATSAT,  SSMSAT  and  WINDSAT  to  the  nature  run,  we  see  that  all  three  analyses 
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are  qualitatively  in  agreement  with  the  nature  run.  Magnitudes  of  the  [v]  are 
all  similar.  However  STATSAT  has  more  small  scale  features  In  the  tropopause 
region  and  SSMSAT  misses  the  strong  northerly  surface  flow  over  the  Southern 
Ocean  at  70° S.  At  the  surface  WINDSAT  and  STATSAT  appear  to  be  roughly 
equivalent . 

Considering  the  amount  of  high  quality  wind  data  available  to  WINDSAT, 
the  small  impr ovement.s  to  the  zonally  averaged  wind  fields  are  disappoint¬ 
ing.  The  lack  of  improvement  in  the  tropical  mean  meridional  circulation  may 
be  caused  by  deficiencies  in  the  assimilation  system.  First,  the  analysis 
wind  increments  must  be  non-divergent :  although  the  height  and  wind  analysis 
decouple  in  the  tropics,  the  wind  analy.sls  still  uses  structure  functions 
derived  from  the  height  structure  function  and  the  assumption  of  geostrophy. 
Secondly,  in  the  NUnHI ,  there  is  no  account  of  the  effects  of  cumulus 
convection. 

The  [T]  has  a  broad  maximum  between  20° S  and  20°N  at  all  levels  througii 
the  tropopause.  Poleward  of  20°,  temperature  decreases.  The  magnitude  of  the 
poteward  temperature  gradient  decreases  with  elevation  up  to  the  200  mb  level 
where  the  gradient  reverses.  Tropopause  height  varies  from  200  mb  in  the 
polar  regions  to  50  mb  at  the  equator.  The  the  nature  run  [T]  agrees  yell 
with  the  GFDL  values  except  that  the  GFDL  tropical  tropopause  is  roughly  at 
the  100  mb  level.  Compared  to  the  nature  run,  we  see  that  STATSAT  is  too  cold 
(by  2  K)  in  the  tropical  PBL  (below  about  900  mb)  and  too  warm  at  the  equator 
at  850  mb  (by  l.A  K) .  This  implies  the  equatorial  region  is  too  stable.  In 
the  S.H.  between  60° S  and  80° S  STATSAT  is  too  cold  below  850  mb  (up  to  -5  K) 
and  there  are  large  positive  errors  over  Antarctica.  Note  that  the  surface  is 
about  700  mb  south  of  80°S.  In  the  Artie  below  850  mb  STATSAT  is  also  too 
warm.  There  are  large  errors  above  the  tropopause;  the  poles  are  warm  and  the 
equatorial  region  is  cold.  WINDSAT  has  errors  similar  to  STATSAT,  however  the 
upper  level  errors  are  substantially  reduced  in  magnitude,  the  equatorial  PBL 
is  even  colder  (by  3  K)  and  the  errors  over  Antarctica  are  reduced.  For  SSM 
the  analyses  are  too  cold  below  850  mb  at  all  latitudes.  In  the  tropics  the 
errors  are  larger  (-5  K)  and  therefore  the  atmosphere  is  even  more  stable. 

The  SSM  analyses  are  very  cold  at  80°N  (250  K) .  Like  WINDSAT  the  errors  above 
the  tropopause  are  smaller  than  in  STATSAT. 
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The  [RH]  in  the  nature  run  (Fig  1 .2h)  is  very  moist  at  the  surface  where 
it  has  a  maxiraun!  of  85%  at  lO'N  and  a  minimum  of  75%  at  30'’N.  [RH]  tends  to 
decrease  from  the  surface  to  500  mb  and  then  increase  again  up  to  300  mb.  In 
addition,  moisture  is  carried  upwards  by  the  ascending  Hadley  circulation  near 
the  equator  and  dry  air  is  brought  downward  by  the  descending  branches  of  the 
Hadley  cell  near  30°N  and  S.  As  a  result,  a  primary  feature  is  the  two 
subtropical  minima  which  occur  in  the  middle  troposphere.  The  northern 
hemisphere  minimum  (about  34%)  occurs  at  500  mb  while  the  southern  minimum 
(31.2%)  lies  near  650  to  700  mb.  This  pattern  is  repeated,  but  with  smaller 
amplitude  by  the  Ferrel  cells. 

[RH]  errors  are  relatively  large  in  all  assimilations.  We  may  contrast 
nature  with  the  cross  sections  for  STATSAT  and  SSMSAT.  In  STATSAT  the  low 
level  averaged  relative  humidity  analysis  below  850  mb  is  consistently  too  dry 
by  10  to  15%.  At  higher  levels  equatorial  relative  humidity  is  too  high  by  5 
to  15%.  And  at  higher  latitudes,  the  upper  level  analysis  is  again  too  dry, 
exceeding  15%  poleward  of  70° .  A  major  result  of  these  analysis  differences 
is  that  the  northern  and  southern  hemisphere  mid- level  minima  are  greatly 
increased  in  magnitude  to  *9.5  and  38.2%,  respectively.  Additionally,  the 
asymmetry  seen  in  nature  with  respect  to  height  is  gone  and  both  features  now 
occur  at  650  mb.  In  SSMSAT  the  averaged  low  level  analysis  is  still  too  dry, 
but  the  difference  is  smaller,  often  less  than  5%.  At  higher  levels  the 
averaged  analysis  is  also  improved  with  differences  less  than  or  equal  to 

5%.  One  exception  is  the  SSMSAT  analysis  at  850  mb  near  30' S  where  the 

difference  with  nature  is  13.8%  compared  with  9.5%  in  STATSAT.  As  in  STATSAT, 
the  averaged  SSMSAT  analysis  does  not  retain  the  asymmetry  in  the  moisture 
field  which  is  seen  in  averaged  nature  data,  although  the  magnitudes  of  the 
minima  are  better  analyzed.  In  short,  as  determined  from  differences  in 
averaged  analyses,  SSMSAT  relative  humidity  analyses  are  closer  to  nature  at 

most  latitudes  and  at  all  vertical  levels.  However,  in  general,  the  polar 

regions  and  boundary  layer  are  too  dry  and  the  mid  latitude  and  tropical 
atmosphere  above  the  PEL  is  too  moist.  Since  moist  air  is  more  buoyant  and 
low  level  moisture  represents  an  energy  source  for  moist  convection,  these 
errors  tend  to  stabilize  the  analyses.  Both  the  temperature  and  humidity 
errors  in  the  tropics  suggest  that  too  much  stabilizing  convection  may  be 
taking  place  in  the  AFGL  model.  The  temperature  errors  might  also  be  due  to 
the  lack  of  a  radiation  parameterization. 
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Relative  humidity  Is  dependent  not  only  upon  the  actual  water  vapor 
amount  but  also  on  the  ambient  air  temperature  which  determines  the  saturation 
value.  If  we  make  the  assumption  that  the  analysis  of  water  vapor  in  its 
conserved  form  (e.g.  mixing  ratio,  specific  humidity)  is  constant  from  nature 
to  the  analysis,  then  analysis  -  nature  differences  in  the  averaged  relative 
humidity  would  tend  to  be  inversely  correlated  with  analysis  -  nature 
differences  in  the  averaged  temperatures. 

For  STATSAT ,  with  the  exception  of  a  broad  area  of  negative  temperature 
differences  over  the  equator  at  500  mb,  we  do  not  find  a  strong  relationship 
between  the  two  variables.  In  fact,  in  some  areas  the  temperature  and 
humidity  differences  seem  to  be  positively  correlated.  In  the  SSMSAT  results 
the  relationship  is  perhaps  stronger  with  a  larger  area  of  negative  tempera¬ 
ture  differences  wiuch  arcs  from  the  surface  in  the  southern  midlatitudes  to 
the  equatorial  midtroposphere  coinciding  with  positive  relative  humidity 
differences  of  order  5  to  10%.  Nevertheless,  the  temperature  differences  tend 
to  be  less  chan  l.S  K  in  most  areas  of  large  rela^’ive  humidity  analysis 
error.  We  cannot  Chen  conclude,  at  least  from  the  zonal- time  averaged 
statistics,  that  most  of  the  errors  in  the  analyzed  humidity  are  due  to 
discrepancies  in  the  temperature  analyses. 

7. A  Calibration  with  OSSE  results 

A  luxury  of  OSSEs  is  the  ability  to  exactly  compute  error  measures .  We 
took  advantage  of  this  in  our  discussion  of  the  analyses.  Analysts  errors  in 
the  real  world  are  not  well  known.  In  fact,  recently  Daley  and  Mayer  (1986) 
presented  analysis  error  of  the  OSSE  experiments  of  Atlas  et  al .  (1985)  (which 
were  discussed  in  the  introduction)  as  surrogates  for  real  analysis  errors. 

We  now  turn  to  an  examination  of  the  forecast  errors.  In  order  to  have  a 
closer  correspondence  with  the  real  world  and  to  simplify  our  calibration 
procedure  we  have  calculated  rms  difference  between  the  forecasts  and  the 
simulated  radiosondes.  We  then  describe  a  procedure  to  calibrate  these 
differences  using  the  NOSAT  -  STATSAT  impact  observed  in  our  previous  OSEs  as 
a  yardstick,  and  present  some  results  of  the  calibration  procedure.  Some 
discussion  of  the  results  of  the  companion  OSSE  experiments.  WINDSAT  is 
included  here  for  comparison. 
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7.4.1  Rms  difference  between  OSSE  forecasts  and  radiosondes 


Rms  differences  between  OSSE  forecasts  and  the  simulated  radiosondes  used 
in  the  data  assimilation  experiments  were  calculated  for  different  regions  and 
for  several  variables  at  each  layer  in  the  atmosphere.  The  variables  examined 
include  geopotential  height,  temperature,  vector  wind,  relative  humidity,  and 
cloud  cover.  Mean  differences  were  also  calculated  and  examined  but  were 
small  relative  to  the  rms  difference  for  nearly  all  variables.  Some  regulari¬ 
ties  observed  in  the  mean  differences  are  described  in  Section  7.4.4.  One 
aspect  of  the  procedure  we  used  is  that  the  forecast  heights  are  anchored  by 
the  verifying  radiosonde  report.  Consequently,  height  errors  described  here 
are  actually  thickness  errors.  As  described  below  we  curve  fit  the  data  to 
determine  impacts  in  terms  of  predictability  time,  l.e.  the  length  of  the 
useful  forecast. 

Due  to  the  variable  density  of  radiosonde  coverage,  global  averages  are 
very  similar  to  N.H.  extratropical  averages  and  S.H.  extratropical  averages 
are  based  on  fairly  small  samples.  Of  course,  these  statistics  are  biased 
towards  land  areas.  Consequently,  differences  between  NOSAT  and  STATSAT  are 
less  pronounced  than  they  might  otherwise  be.  As  a  result  of  these  factors  as 
well  as  the  radiosonde  errors  themselves,  we  expect  only  qualitative  ^reement 
between  impacts  described  here  and  impacts  measured  in  Section  4.2  by 
comparing  forecasts  grid  point  by  grid  point  to  the  nature  run. 

Examples  of  the  growth  of  forecast  errors  as  evidenced  by  the  rms  dif¬ 
ference  are  displayed  in  Figs.  7.25  through  7.27.  In  each  figure  there  are 
three  panels,  showing  results  averaged  over  3  forecasts  for  the  N.H. 
extratropics,  tropics  and  S.H.  extratropics.  Here  the  tropics  are  taken  to 
run  from  30* S  to  30*N.  Rms  difference  curves  for  smaller  regions  which  were 
studied,  generally  behaved  as  described  here  for  the  larger  parent  regions. 

In  the  S.H.,  the  rms  difference  curves  sometimes  exhibit  a  sawtooth 
pattern  due  to  sampling  problems;  there  are  usually  about  60  RAOBs  at  00  GMT 
and  only  about  40  at  12  GMT  in  the  S.H.  Most  of  the  non- reporting  RAOBs  are 
in  the  Australian  sector. 

Considering,  first  the  rms  difference  for  500  mb  geopotential  (Fig.  7.25) 
we  see  that  impacts  in  the  N.H.  are  relatively  small.  WINDSAT  lags  STATSAT  by 
approximately  7  h.  Midway  between  these  two  lie  SSMSAT  and  NOSAT.  SSM+TOVS 
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IS  equivalent  to  STATSAT  In  the  tropics,  the  differences  are  all  near  20  m 
and  the  growth  rate  is  very  small.  Apparently,  the  tropical  analysis  error;, 
are  so  large  chat  the  error  growth  is  near  saturation.  At  least  during  the 
initial  part  of  the  forecast  WINDSAT  has  a  slight  edge,  which  because  of  the 
slow  growth  of  differences  corresponds  to  9  or  10  hours  of  predictability 
time.  Impacts  in  the  S.H.  are  very  large.  WINDSAT  is  36  h  better  than 
STATSAT,  which  is  in  turn  more  than  36  h  better  than  the  NOSAT  forecasts.  The 
three  sounder  based  systems  are  roughly  equivalent  with  SSMSAT  and  SSM+TOVS 
tetter  chan  STATSAT  by  12  and  8  hours  respec  t ive ly . 

Next,  we  consider  the  rms  vector  wind  differences  at  200  mb 
(Fig.  7,26),  Clearly,  WINDSAT  always  yields  a  big  improvement.  Compared  to 
ST.MSAT,  WINDSAT  provides  1,  2  and  2.75  day  improvements  in  forecast  skill  in 
ihc  N.H.,  tropics  and  S.H,  respectively.  NOSAT  is  particularly  poor  ovei  ilic 
S .  ii .  and  the  tropics.  It  is  quite  good  over  Europe,  where  it  outpcrfoj'iii.s  hotli 
SSN  +  TCVo  and  STATSAT.  Of  the  three  sounder  syscem.s,  SSMSAT  is  generally 
better,  improving  predictability  by  at  least  a  day  in  the  S.H.  and  tropic;;. 

The  rms  difference  for  relative  humidity  at  850  mb  are  shown  in 
tig.  7.27.  Impacts  in  terms  of  forecast  time  are  all  relatively  small  in  tiie 
N.H.,  generally  in  the  range  3  to  12  h.  In  the  extratropics,  WINDSAT, is  the 
best.  This  might  have  been  anticipated  since  the  WINDSAT  analyses  of  mass  and 
wind  are  superior  in  the  extratropics  and  since  relative  humidity  is  so 
strongly  influenced  by  the  large  scale  synoptic  systems  which  are  better 
forecast  by  WINDSAT.  In  the  tropics,  the  SSMSAT  moisture  forecasts  arc  the 
best.  Overall  the  ranking  is  WINDSAT,  SSMSAT,  NOSAT,  SSM+TOVS  and  STATSAT. 

It  appears  that  using  TOVS  degrades  the  moisture  analysis.  Since  TOVS 
relative  humidity  retrievals  were  not  used,  this  poses  a  conundrum.  TOVS  data 
affect  the  model  specific  humidity  Indirectly  because  the  analyzed  'unriablcs 
are  temperature  and  relative  humidity.  That  is,  in  an  area  with  only  TOVS 
data  the  updated  temperature  field  is  combined  with  the  unaltered  relative 
humidity  field  to  update  the  model  specific  hujnidity. 

We  also  calculated  rms  difference  in  cloud  cover  layer  by  layer.  Cloud 
cover  was  diagnosed  from  relative  humidity  using  the  inverse  Tibaldi  scliemc 
(as  described  in  Section  6.2).  Invariably,  the  corresponding  relative 
humidity  and  cloud  cover  plots  look  very  similar.  Fig.  7.28  shows  the  glob.il 
rtTi.s  difference  for  850  mb  relative  humidity  and  cloud  cover.  Except  ior  tliu 
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fact  that  the  cloud  cover  errors  are  larger  In  magnitude,  the  curves  in  the 
two  panels  of  Fig.  7.28  are  nearly  the  same.  For  this  reason  we  have  not 
included  any  other  cloud  cover  rms  difference  statistics  in  this  report 

We  examined  the  growth  of  rms  difference  for  levels  others  than  those 
described  here.  The  results  shown  here  are  generally  representative.  In  the 
next  section  tables  summarizing  all  the  levels  are  given.  One  notable 
exception  to  the  general  trend  discussed  so  far,  is  the  positive  impact  of 
SSM+TOVS  on  the  upper  level  heights.  For  example,  Fig.  7.29  shows  the  N.H. 

50  mb  height  rms  difference  evolution.  Here  the  ordering  is  SSM+TOVS, 

STATSAT,  WINDSAT  and  NOSAT  with  an  Increment  of  about  6  hour  s  between  each 
pair;  SSMSAT  and  STATSAT  are  roughly  equivalent. 

7.4.2  Calibration  procedure 

It  has  been  observed  that  OSSE  forecasts  are  too  good  because  any  two 
models,  such  as  the  model  used  in  the  experiments  and  the  model  used  Co 
generate  nature,  are  more  alike  than  any  model  and  the  real  atmosphere. 
Consequently,  it  is  unwise  to  naively  carry  over  the  forecast  impacts  observed 
in  OSSEs  to  the  real  world.  For  example,  at  short  forecast  times,  OSSE 
forecasts  tend  to  be  so  good  that  there  is  little  room  for  improvement,;  adding 
a  new  observing  system  might  then  have  a  smaller  impact  than  in  the  real 
world.  On  the  other  hand,  at  longer  forecast  times,  real  data  forecasts  will 
be  so  bad  that  a  new  observing  system  will  have  no  Impact  while  the 
corresponding  OSSE  Impact  may  be  significant. 

For  these  reasons  it  is  desirable  to  calibrate  the  OSSE  results. 

However,  for  the  present  experiments  we  find  that  the  OSSE  impacts  are  fairly 
similar  to  the  OSE  impacts  and  the  calibration  procedure  does  not  greatly 
alter  the  conclusions  one  might  draw  from  the  OSSE  results  directly.  To 
minimize  practical  and  interpretive  difficulties  we  use  only  very  simple 
calibration  approaches.  Our  principal  assumption  is  that  the  OSE  impa..L  of 
adding  or  removing  an  observing  system  is  proportional  to  the  corresponding 
OSSE  impact.  In  our  calibrations  we  always  take  STATSAT  to  be  our  standard. 

We  use  the  NOSAT  -  STATSAT  differences  to  determine  the  constant  of 
proportionality . 

Impact  may  be  measured  in  many  ways.  Useful  impact  measures  should 
account  for  differences  in  variability  from  season  to  season  and  from  region 
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Lo  i-' .  example  an  impact  of  10  geopotential  meters  is  meaningless 

without  t!  o  context  of  place,  season  and  vertical  level.  Usually,  for  the 
pvirpc  .e  of  comparison,  it  is  reasonable  to  scale  the  squared  errors  by  ttu  ii 
respective  climate  variance. 

It  i.s  al.yo  often  helpful  to  define  impact  In  terms  of  the  change  i  r. 
prc'-lictahi  1  itv  ti;;ie  For  example,  one  might  define  the  predictability  time  ,c, 
t'ne  tiite  at.  which  the  mean  squared  forecast  error  reaches  the  climate  v.ii' i  .irce 
level.  Foreca.sts  with  errors  this  large  are  normally  worthless.  A  posit  i 
impact  it',  predictability  time  would  then  indicate  the  adciitional  time  that  :  lie 
tuttcast  remains  useful. 

Measuring  imp.ict  in  terms  of  p redi c tabi 1 i ty  time  is  especially  useful 
when  cerbined  with  a  common  idealization  of  the  growth  of  forecast  error.  .A 
nuithe;  ot  s  i  ii.p  1  e  pa  rattle  to  r  iza  c  ions  of  the  growth  of  error  have  been  advaiu  t  d 
(o.g.  Ltilclier  ai,d  Kalntiv,  1987  and  references  therein).  Remarkably  good  1  i  i  :> 
to  en.semble  averaged  forecast  error  growth  curves  have  been  obtained,  by 
fitting  relatively  simple  autonomous  constant  coefficietnt  ordinary 
di f feretit ial  equations.  These  coefficients  describe  the  growth  of  small 
.■rrors,  the  saturation  of  large  errors  at  the  climate  variance  level  and  the 
source  of  errors  due  to  modeling  deficiencies.  Since  these  constantrs  should 
be  the  same  for  a  set  of  experiments,  e.g.  for  all  our  February  OSEs,  all 
fitted  error  growth  curves  for  the  experiments  should  be  the  same  except  for  a 
translation  with  re.specr  to  the  time  axis.  This  shift  is  the  impact  in  term 
of  predictability  time. 

A  direct  reading  of  the  predictability  times  from  the  rms  difference 
curves  proved  difficult  because  our  sample  is  rather  small.  We  could  fit  tlu 
parameterization  of  Dalcher  and  Kalnay  (1987).  Instead,  we  took  advantage  of 
the  observation  that  our  rms  difference  curves  grow  nearly  linearly,  at  least 
during  the  forecast  period  from  12  to  98  hours,  to  fit  these  data  with  a 
series  of  straight  lines  having  a  common  slope.  In  the  N.H.  these  fits  were 
very  good.  They  are  less  reliable  in  the  S.H.  and  tropics,  presumably  because 
the  number  of  radiosondes  in  these  regions  ’s  small.  Our  results  for  heigtu, 
wind  and  relative  humidity  are  shovm  In  Tables  7.1  through  7.3.  In  the 
tafjles,  each  item  in  the  columns  labeled  R**2  is  the  fraction  of  variance 
explained  by  the  fitting  procedure  or  equivalently  the  square  of  the  coriel.i- 
tion  coefficient.  The  other  columns  in  the  tables  display  the  predic i ah i 1 i t y 
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time  impacts  of  the  various  experiments,  all  relative  to  STATSAT.  Some  of 
these  values  have  been  quoted  earlier.  These  impacts  are  the  difference  in 
the  X  (or  time)  intercepts  of  the  fitted  lines.  (These  intercepts  are 
proportional  to  the  y  (or  rms  difference)  intercepts  with  proportionality 
factor  equal  to  the  common  slope.) 

The  NOSAT  impact,  i.e.  the  difference  between  NOSAT  and  STATSAT,  in  the 
OSE  experiments  were  then  used  to  calibrate  the  OSSE  results  according  to 

(Expected  OSE  impact)  -  (NOSAT  OSE  impact)*(OSSE  impact)/ 

(NOSAT  OSSE  impact) 

This  provided  us  with  calibrated  intercepts  which  we  combined  with  the 
observed  OSE  error  growth  rate  (the  common  slope)  to  create  calibrated  rms 
difference  curves  for  WINDSAT,  SSM  and  SSM+TOVS  experiments.  These  are 
displayed  along  with  the  observed  OSE  results  in  the  figures.  The  horizontal 
distances  between  the  various  curves  are  the  predictability  time  impacts. 

7.4.3  Calibration  results 

The  calibration  procedure  described  above  allows  us  to  translate  our  OSSE 
results  into  anticipated  real  world  impacts  in  a  quantitative  manner.  Several 
examples  are  provided  in  Figs.  7.30  through  7.32.  Complete  details  are 
provided  by  Tables  7.1  through  7.3  and  the  above  equation  for  the  expected 
impact . 

In  Fig.  7.30,  the  S.H.  500  mb  height  rms  difference  for  the  STATSAT  and 
NOSAT  OSE  experiments  are  plotted.  The  calibrated  OSSE  results  for  SSM, 
SSM+TOVS  and  WINDSAT  are  also  plotted.  These  are  the  three  straight  lines 
between  12  and  48  h  on  the  plots.  They  are  plotted  only  for  this  period  since 
it  is  only  this  period  which  was  used  in  the  curve  fitting.  Note  that  the 
February  (a)  and  June  (b)  OSEs  provide  two  independent  calibrations.  In  both 
cases  the  dramatic  improvement  seen  in  the  OSSEs  for  WINDSAT  is  expected  to 
carry  over  in  actuality.  A  36  h  improvement  in  forecast  skill  relative  to 
STATSAT  is  anticipated  and  the  expected  impact  of  the  SSM  data  is  12  h.  In 
the  N.H.  (not  shown)  STATSAT  and  NOSAT  OSSE  results  are  nearly  equivalent,  so 
impacts  expected  from  advanced  observing  systems  cannot  be  calibrated. 


96 


The  tropical  200  mb  vector  wind  rms  differences  are  shown  for  February 
(a)  and  June  (b)  OSSEs  (Fig.  7.31).  Because  the  growth  rates  for  the 
differences  are  so  small,  it  is  difficult  to  judge  the  impacts  by  eye.  In  this 
case,  in  terms  of  predictability  time  there  are  significant  impacts.  WINDSAT 
has  a  24  h  advantage  over  SSMSAT,  a  48  h  advantage  over  SSM+TOVS  and  a  60  h 
advantage  over  STATSAT.  Further  SSM+TOVS  has  approximately  a  12  h  advantage 
over  STATSAT  which  in  turn  has  a  12  h  advantage  over  NOSAT.  However,  not  much 
weight  should  be  given  these  results  since  in  the  tropics  and  S.H.,  the 
califc>ration  procedure  is  not  very  reliable. 

For  the  N.H.  850  nib  relative  humidity  or  cloud  cover  (Fig.  7.32),  the 
actual  impacts  are  all  expected  to  be  rather  small.  Note  again  the  close 
relationship  between  relative  humidity  and  cloud  cover  statistics. 

TVie  calibration  procedure  does  have  some  uncertainties  and  drawbacks. 

The  main  drawback  is  the  assumption  of  a  linear  relationship  between  impact  in 
the  OSSEs  and  in  the  OSEs.  Of  course  the  data  assimilation  system  and  nature 
are  highly  nonlinear.  In  the  current  experiments,  the  uncertainties  are 
mostly  due  to  the  small  sample  size,  especially  in  the  S.H.  where  the  number 
of  radiosondes  used  in  the  verification  is  small.  For  example,  in  some  cases 
the  sense  of  impact  between  NOSAT  and  STATSAT  is  reversed  in  OSSEs  ahd^OSEs. 

In  these  cases  the  calibration  produces  nonsensical  results.  This  occurred 
when  calibrating  the  S.H.  850  mb  relative  humidity  rms  difference.  In  other 
cases  the  OSE  impact  between  NOSAT  and  STATSAT  is  quite  small.  This  implies 
negligible  calibrated  impact  for  any  change  to  the  data  assimilation  system. 
This  occurred  when  calibrating  the  N.H.  500  mb  height  rms  difference. 

7.4,4  Forecast  biases 

In  general  the  biases  during  the  forecast  are  small  compared  to  the  rms 
difference.  However  in  many  cases  the  biases  grow  very  steadily  with  time 
indicating  that  the  AFGL  model  is  warming  and  drying  relative  to  the  ECMWF 
nature . 

We  examined  ''be  biases  by  fitting  straight  lines  with  a  common  slope  as 
described  in  Section  7.4.2.  However  in  this  case  all  data  from  12  to  96  hours 
was  used.  The  common  slope  obtained  from  the  fit  is  then  the  rate  at  which 
the  biases  grow.  Some  of  these  are  displayed  in  Table  7.4  for  the  N.H. 

OSSEs.  Generally  speaking  the  forecasts  are  warming  relative  to  nature  by  one 
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third  to  one  degree  per  day.  The  height  biases  reflect  these  temperature 
biases.  These  results  are  consistent  with  the  warming  seen  during  the  spinup 
forecast  (Section  4.1)  Also  the  forecasts  are  drying  in  the  lower  atmosphere 
by  1  to  1.5  percentage  points  of  relative  humidity  per  day.  For  temperature 
and  height  there  are  many  cases  when  the  fraction  of  variance  explained  by  the 
fit  is  greater  than  .99,  indicating  that  Che  bias  grows  very  linearly.  For 
example  Fig.  7.33  shows  the  evolution  of  bias  for  the  N.H.  500  mb  height  in 
the  OSSEs .  Differences  between  the  experiments  are  not  significant.  Results 
for  the  tropics  and  S.H.  are  not  as  regular  and  clear  cut  presumably  because 
of  sampling  variability.  This  also  applies  to  the  OSEs ,  although  there  is 
some  evidence  of  the  forecasts  warming  during  the  February  OSEs  in  the  mid 
troposphere.  In  the  N.H.  the  dry  relative  humidity  biases  are  substantial  at 
1000  mb.  In  this  case  the  bias  at  the  start  of  Che  forecast  is  already  -17  to 
-18  percent.  At  other  levels  the  initial  dry  bias  is  only  of  order  5%. 
Typically,  the  analyses  are  dry  by  5%  and  the  forecasts  continue  to  dry  out  by 
1%  to  2%  per  day  for  the  first  two  days  of  the  forecast.  For  example,  the 
evolution  of  the  biases  of  the  850  mb  relative  humidity  forecasts  are  shown  in 
Fig.  7.34.  Again,  differences  in  the  experiments  are  not  significant.  The 
trends  are  not  so  clear  cut  in  the  tropics  and  extratropics  and  the  OSEs.  In 
the  tropics,  in  the  OSSEs,  the  atmosphere  tends  to  moisten  during  the  *. 
forecast.  In  this  case  the  boundary  layer  is  analyzed  dry  but  the  mid 
troposphere  is  slightly  moist. 
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Table  ^.1  Predictability  Impacts  (days)  for  height  (a),  wind  (b)  and  relative 
humidity  (c)  for  the  N.H.  extratropics.  Impacts  for  cases  with 
R**2  less  than  0.25  (i.e.  for  correlations  less  than  0.5)  are  not 
shown . 


A.  Height 

SSSSSSSSSSSSSSSXSXSSXSSSSCSXSSSSSSKSSSSSSXeCSCSSXSSSBSSSSKSSSSSSnMSKKSSSBSSSBSSSBSSSSnsSSCSSSXS 


Noventier  OSSE 


February  OSE  Jirte  OSE 


I  P{1) 


R**2 

NOSAT 

SSH/T 

SSM+TOVS 

UtNDSAT 

R**2 

NOSAT 

R**2 

NOSAT 

1 

1000 

.000 

- 

- 

- 

.603 

.266 

.219 

- 

2 

850 

.822 

.183 

-.143 

.022 

.601 

.975 

-.287 

,920 

-.319 

3 

700 

.991 

.168 

.118 

,011 

.488 

.996 

-.048 

.989 

-.032 

i 

500 

.999 

.U1 

.1U 

-.008 

.307 

.999 

.023 

.985 

-.003 

5 

^00 

.999 

.098 

.113 

-.017 

.255 

.998 

.022 

.960 

o 

o 

6 

300 

.998 

.038 

.064 

-.010 

.186 

.998 

.000 

.978 

-.070 

7 

250 

.996 

.019 

.050 

-.013 

.167 

.998 

-.019 

.958 

-.115 

8 

200 

.996 

-.011 

.045 

.012 

.103 

.998 

-.051 

.874 

-.145 

9 

150 

.994 

-.013 

.060 

.029 

.120 

.997 

-.104 

.778 

-.175 

10 

100 

.996 

-.100 

.049 

.065 

-.112 

.993 

-.121 

.449 

-.444 

11 

70 

.999 

-.412 

.060 

.214 

-.294 

.990 

-.127 

.111 

- 

12 

50 

.998 

-.699 

.064 

.307 

-.368 

.972 

-.220 

.434 

-.585 

B.  Vector  wind 

SBSSSSSSSSSSBSSSSSSSSSSSSBBSSSBSSSSSSSSSSSZBSSSSSSSSSSSBSSSBSSSSSBSSCSSBSSSSSSBSSSrsSSrXSSXBSSirBS 


I  P(l) 


Novenber  OSSE 


February  OSE  June  OSE 


R»»2 

NOSAT 

SSN/T 

SSM+TOVS 

Wl  NOSAT 

R**2 

NOSAT 

R**2 

NOSAT  . 

1 

1000 

.969 

.414 

.191 

-.094 

.928 

.909 

.161 

.910 

.451 

2 

850 

.963 

.474 

.207 

-.210 

.915 

.987 

.122 

.977 

-.144 

3 

700 

.981 

.458 

.149 

-.314 

.973 

.990 

.040 

.932 

-.184 

4 

500 

.976 

.410 

.179 

-.176 

.866 

.996 

.015 

.948 

-.160 

5 

400 

.985 

.373 

.182 

-.145 

.846 

.995 

.009 

.984 

-.152 

6 

300 

.980 

.233 

.134 

-.154 

.806 

.977 

-.056 

.990 

-.103 

7 

250 

.970 

.214 

.147 

-.163 

.941 

.989 

-.060 

.987 

-.135 

8 

200 

.948 

.126 

.214 

-.142 

1.079 

.979 

-.171 

.990 

-.235 

9 

150 

.860 

.001 

.318 

-.095 

1.301 

.965 

-.292 

.988 

-.250 

10 

100 

.892 

.232 

.489 

.127 

1.244 

.929 

-.175 

.984 

-.108 

11 

70 

.884 

.265 

.572 

.331 

1.222 

.893 

-.035 

.964 

.010 

12 

50 

.929 

.245 

.506 

.399 

1.001 

.938 

-.076 

.970 

.157 

C.  Relative  hunidity 

SSSSSSSSSSSSBSXSSXBSasSSSSSSSSCSSSSSXSSSSeSSSSSSSSSSSBSeSSCSXBSKBSBBBSCeXXSSSSSSCSSSSSSSSSSSSSSSS 


Novenber  OSSE  February  OSE  June  OSE 


1 

P(l) 

R**2 

NOSAT 

SSM/T 

SSM+TOVS 

Wl  NOSAT 

R**2 

NOSAT 

R**2 

NOSAT 

1 

1000 

.781 

.213 

.197 

-.013 

.490 

.822 

-.154 

.107 

- 

2 

850 

.959 

.320 

.388 

.168 

.543 

.935 

.100 

.964 

.233 

3 

700 

.954 

.365 

.330 

-.029 

.719 

.856 

.136 

.945 

.120 

4 

500 

.906 

.475 

.392 

.242 

1.121 

.973 

.400 

.934 

.262 

5 

400 

.921 

.161 

.116 

.014 

.763 

.990 

.114 

.886 

-.019 

6 

300 

.853 

.188 

.146 

.066 

.906 

.861 

.355 

.846 

-.086 

99 


Table  7.2  Predictability  impacts  (days)  for  height  (a),  wind  (b)  and  relative 
humidity  (c)  for  the  tropics.  Impacts  for  cases  with  R**2  less 
than  0.25  (i.e.  for  correlations  less  than  0.5)  are  not  shown. 


A.  Height 

1  P(t) 

BSSSSBSB888SSXSSCSI 

N(r 

««^r  OSSE 

February  OSE  June  OSE 

R**2 

NOSAT 

SSM/T  SStHTOVS 

WINDSAT 

R**2 

NOSAT 

R**2 

NOSAT 

1 

1000 

.000 

. 

. 

- 

. 

.158 

. 

.312 

.181 

2 

850 

.131 

- 

- 

- 

- 

.126 

- 

.643 

-.466 

3 

700 

.652 

-.490 

-.033 

.176 

.U5 

.484 

.395 

.779 

-.641 

4 

500 

.551 

-.224 

-.145 

-.623 

.610 

.028 

- 

.241 

- 

5 

400 

.569 

-.348 

-.302 

-1.007 

.398 

.695 

.066 

.672 

-.804 

6 

300 

.381 

-.687 

-.364 

-.917 

.414 

.633 

.658 

.545 

-.800 

7 

250 

.447 

-.499 

.116 

-.050 

.465 

.000 

.880 

-.921 

fi 

200 

.411 

-.147 

.731 

.665 

.940 

.001 

.881 

-.517 

9 

150 

.191 

- 

- 

- 

- 

.041 

.962 

-.348 

10 

100 

.171 

- 

- 

- 

- 

.130 

.556 

-.354 

11 

70 

.009 

- 

- 

- 

- 

.094 

.005 

- 

12 

50 

.005 

- 

- 

- 

- 

.002 

.000 

- 

B.  Vector  wind 


ssssssxxssssssssssssasesssssssss 


I  P<1) 


Novealier  OSSE 


=SSX=SSBCX=asSSB=S£SSSCSSrXSCS&S=SSS 

February  OSE  June  OSE 


R**2 

NOSAT 

SSH/T 

SSM+TOVS 

WINDSAT 

R**2 

NOSAT 

R**2 

NOSAT 

1 

1000 

.450 

.802 

.815 

.458 

1.389 

.047 

• 

.641 

.323 

2 

850 

.846 

.416 

.082 

.029 

.846 

.813 

-.246 

.791 

-.304 

3 

700 

.881 

.771 

.765 

.479 

1.307 

.368 

.968 

.883 

.165 

4 

500 

.550 

.882 

.956 

.267 

1.682 

.402 

.023 

.770 

-.321 

5 

400 

.565 

.861 

1.150 

.664 

2.000 

.334 

.472 

.775 

-.794 

6 

300 

.494 

.261 

1.762 

-.128 

3.908 

.806 

-.433 

.537 

-.575 

7 

250 

.398 

.418 

1.454 

.693 

3.479 

.788 

-.212 

.638 

-.089 

8 

200 

.664 

-.462 

1.317 

.482 

2.057 

.666 

-.649 

.345 

-.449 

9 

150 

.739 

.168 

.932 

.350 

2.624 

.328 

1.527 

.653 

-.433 

10 

100 

.033 

- 

- 

- 

- 

.577 

-.536 

.902 

-.158 

11 

70 

.873 

.355 

.854 

.636 

1.719 

.015 

- 

.736 

-.267 

12 

50 

.698 

1.160 

2.006 

1.534 

3.665 

.049 

- 

.018 

- 

C.  Relative  hwnidity 

SSSSSSSSSSSSSCSSSXXSBXSSSSS8SS«BSSSBSCX8SXS8BSS3SSSSSXSSSSSSSSSSSseSSSSSSS8 


Noveaber  OSSE 


February  OSE 


B8SS8BBS8S8SSSSSBS 

June  OSE 


1 

P<1) 

R»»2 

NOSAT 

SSM/T  SSM+TOVS 

WINDSAT 

R*»2 

NOSAT 

R»*2 

NOSAT 

1 

1000 

.114 

- 

- 

- 

.256 

-1.478 

.522 

.319 

2 

850 

.424 

1.520 

2.104 

1.775 

1.671 

.521 

-.121 

.944 

.071 

3 

700 

.433 

.601 

.349 

.460 

1.018 

.830 

.210 

.816 

-.285 

4 

500 

.795 

-.237 

.064 

.086 

.149 

.000 

- 

.899 

-.098 

5 

400 

.624 

-.135 

-.192 

-.243 

-.059 

.545 

1.828 

.393 

.016 

6 

300 

.157 

• 

- 

• 

- 

.467 

-.682 

.480 

-.275 

100 


Table  7.3 


Predictabili  ty  Imacts  (days)  for  height  (a)  ,  wind  (b)  and  relative 
humidity  (c)  for  the  S.H.  extratropics.  Impacts  for  cases  with 
R**2  less  than  2.5  (l.e.  for  correlation  less  than  0.5)  are  not 
shown . 


A.  Height 

S==SSS  =  ==S  =  5  =  S  =  SX  =  S  =  =S=£KSSSSaKU«K«CUSXKKKX«XMK«SXK«SKnX«S3BES«CESKXS=rBCSXSBSSSBS«nSnXSS3B 

Novonber  OSSE  February  OSE  June  OSE 


I  P(l) 


R»*2 

NOSAT 

SSH/T 

SSH+TOVS 

UINOSAT 

R**2 

NOSAT 

R**2 

NOSAT 

1 

1000 

.000 

. 

. 

. 

. 

.062 

. 

.024 

- 

2 

850 

.486 

-.082 

2.098 

2.082 

1.954 

.000 

- 

.308 

-.363 

3 

700 

.795 

-1.076 

.469 

.471 

1.357 

.582 

-1.396 

.318 

-1.958 

4 

500 

.671 

-1.597 

.498 

.285 

1.686 

.678 

-1.231 

.688 

-1.372 

5 

400 

.668 

-1.634 

.524 

.323 

1.491 

.469 

-1.625 

.777 

-1.783 

6 

300 

.734 

-1.709 

.472 

.350 

1.106 

.399 

-1.592 

.751 

-2.208 

7 

250 

.742 

-1.659 

.528 

.448 

.949 

.570 

-1.423 

.711 

-3.210 

8 

200 

.747 

-1.605 

.619 

.502 

.784 

.403 

-1.765 

.598 

-4.398 

9 

150 

.717 

-1.768 

.626 

.518 

.798 

.369 

-1.745 

.689 

-4.280 

10 

100 

.647 

-2.243 

.318 

.267 

.193 

.483 

-1.111 

.560 

-5.290 

11 

70 

.577 

-3.387 

-.004 

.841 

-1.169 

.802 

-.787 

.013 

- 

12 

50 

.610 

-4.762 

-.845 

1.256 

-2.376 

.764 

-.337 

.184 

- 

B.  Vector  wind 

B=SX=SSSS;SSSXSSSSSSSS  =  =  X  =  =  SSUBBaSXBmBBBBKBBXSSSBBeBBSSCBS«SSSSBBSSC=SES=SS£S  =  S=:  =  B=  ===£BSSS3S=KBS 


November  OSSE 


February  OSE  June  OSE 


I  PCD 


R*»2 

NOSAT 

SSN/T 

SSM+TOVS 

UtNDSAT 

R**2 

NOSAT 

R**2 

NOSAT 

1 

1000 

.145 
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• 

• 

.057 

. 

.304 

-1.757 

2 

850 

.716 

-.518 

1.674 

1.119 

3.372 

.795 

-1.274 

.688 

-2.197 

3 

700 

.696 

-1.133 

1.117 

.406 

2.559 

.265 

-2.042 

.700 

-2.487 

4 

500 

.651 

-1.435 

.644 

.172 

2.306 

.672 

-.651 

.320 

-2.112 

5 

400 

.601 

-1.615 

.536 

.U5 

2.571 

.573 

-.714 

.454 

-1.461 

6 

300 

.852 

-1.547 

.333 

.081 

2.055 

.242 

- 

.154 

7 

250 

.819 

-1.359 

.135 

.043 

1.516 

.244 

- 

.024 

8 

200 

.808 

-1.562 

1.011 

.539 

2.749 

.777 

-1.113 

.004 

9 

150 

.520 

-1.716 

1.448 

.390 

3.460 

.117 

- 

.016 

10 

100 

.051 

- 

- 

- 

- 

.016 

- 

.124 

11 

70 

.265 

-1.454 

1.621 

.503 

4.040 

.059 

- 

.575 

-3.541 

12 

SO 

.169 

- 

- 

- 

- 

.290 

-.335 

.418 

-1.083 

C.  Relative  huxidity 

«ss8asssgrsgsrgsssssgssgggg»»g«»»«m«»»««««'Bwr— rg«g«gCTc«gg»Kc»Kc»s«»gKgarc«sgssacsgK8gsK«g«wBg»K«« 


Noveatier  OSSE 


February  OSE  June  OSE 


1 

P(l) 

R»»2 

NOSAT 

SSH/T  SSM+TOVS 

UINDSAT 

R*»2 

NOSAT 

R**2 

NOSAT 

1 

1000 

.133 

- 

. 

- 

- 

.005 

- 

.705 

-.075 

2 

850 

.432 

.062 

1.102 

.510 

2.193 

.364 

-.969 

.140 

- 

3 

700 

.741 

.398 

.832 

.314 

1.263 

.188 

- 

.067 

- 

4 

500 

.322 

-.026 

2.434 

.407 

2.980 

.045 

- 

.662 

.225 

5 

400 

.427 

.014 

1.667 

.293 

2.778 

.493 

-.037 

.005 

- 

6 

300 

.168 

• 

• 

- 

- 

.158 

- 

.705 

.664 

101 


Table  7.4  Growth  rate  of  forecast  bias  for  height  (m/day) ,  temperature 

(K/day)  and  relative  humidity  (%/day)  for  the  N.H.  extratropics. 
Impacts  for  cases  with  R**2  less  than  0.64  (i.e.  for  correaltlons 
less  than  0.8)  are  not  shown. 


1 

P(l) 

R**2 

Z 

R**2 

T 

R**2 

RH 

1 

1000 

.000 

* 

.175 

. 

.673 

-1.467 

2 

8S0 

.561 

- 

.879 

.443 

.753 

-  .948 

3 

700 

.963 

4.232 

.983 

.718 

.848 

-1.074 

4 

500 

.989 

12.625 

.991 

.937 

.284 

- 

5 

400 

.992 

18.747 

.986 

.903 

.014 

- 

6 

300 

.992 

25.569 

.986 

.689 

.057 

- 

7 

250 

.993 

28.995 

.985 

.544 

- 

- 

8 

200 

.994 

31.843 

.979 

.426 

- 

- 

9 

150 

.995 

35.529 

.955 

.371 

- 

- 

10 

100 

.996 

39.334 

.970 

.375 

- 

- 

11 

70 

.998 

43.682 

.964 

.406 

- 

- 

12 

50 

.997 

47.277 

.973 

.415 

- 

- 
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Figure  7.1  Northern  hemisphere  height  analyses  and  analysis  errors  at  00  GMT 
23  November.  The  contour  Interval  in  the  500  mb  analyses  is  80  m 
at  1000  mb  it  is  40  m.  The  contour  interval  of  the  error  maps  is 
40  m  at  both  1000  and  500  mb.  Shown  are  500  mb  nature  data  (a). 
1000  mb  nature  data  (b) . 


Figure  7 . 


2  Format  as  in  Fig.  7.1.  STATSAT  analysis  (a)  and  analysis  error 
(b)  at  500  mb,  and  at  1000  mb  (c.d). 
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Figure  7.3  Format  as  in  Fig.  7.1.  SSMSAT  analysis  (a)  and  analysis  error  (b) 
at  500  mb,  and  at  1000  mb  {c,d).  Raob  height  residuals  at  o-0.5 
used  in  the  analyses  are  marked  by  crosses  in  panel  b. 
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Figure  7. 


3  (continued) . 
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Figure  7 . 


4  Format  as  In  Fig.  7.1.  SSM+TOVS  analysis  (a)  and  analysis  error 
(b)  at  500  mb,  and  at  1000  mb  (c,d).  Raob  height  residuals  at 
O-0.5  used  in  the  analyses  are  marked  by  crosses  in  panel  b. 
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Figure  7.5  Northern  hemisphere  height  analyses  and  analysis  errors  at  00  GMl 
25  November.  The  contour  interval  in  the  500  mb  analyses  is  80  m 
at  1000  mb  it  is  40  m.  The  contour  interval  of  the  error  maps  is 
40  m  at  both  1000  and  500  mb.  Shovm  are  500  mb  nature  data  (a'. 
1000  mb  nature  data  (b) . 


rora 


Figure  7.8  (continued). 
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Figure  7.12  Format  as  In  Fig.  7.11.  STATSAT  analysis  (a)  and  analysis  error 
(b)  at  500  mb,  and  at  1000  mb  (c.d). 
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urriTUX  <  oeates  i  udiiu*  <  otcmcs  i  uTinae  i  acucci 


L0C(njK  I  (ccRcn  esf  i 


Figure  7.17  Inferred  low  cloud  cover  from  relative  humidity  valid 
00  GMT  25  November.  (a)  Nature,  (b)  STATSAT  analysis, 

(c)  STATSAT  -  Nature.  Contour  interval  is  25%  with  zero  line 
omitted.  In  difference  field,  negative  values  are  dashed. 
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IHTnuOC  C  MCRteS 


Figure  7.18  Inferred  low  cloud  cover  from  relative  humidity  valid  00  GMT 

25  November.  (a)  SSMSAT  analysis,  (b)  SSMSAT  -  Nature.  Contour 
Interval  is  25%  with  zero  line  omitted.  In  difference  field, 
negative  values  are  dashed. 
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GUf.  iciCHr  mx  <  n  i 


TIK  I  JULim  ORT  I 


Figure  7.19  Rms  height  errors  at  500  mb.  (a)  NOSAT,  (b)  STATSAT,  (c)  SSMSAT, 
(d)  SSM+TOVS.  Analysis  errors  are  shovm  in  solid  curves, 
forecast  errors  in  dashed  curves.  Julian  day  322  corresponds  to 
00  GMT  18  November. 
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MiiiiBrt*  I  1  mLiiHMi 


Figure  7.21  (continued). 


LSTinae  i  acEKCs  i  mnnaE  i  9eckcs  ) 


Figure  7.22  Zonal  time  averaged  v  component  of  wind.  (a)  Nature, 

(b)  STATSAT,  (c)  STATSAT  -  Nature,  (a)  WINDSAT,  (e)  WINDSAT  - 
Nature,  (f)  SSMSAT,  (g)  SSMSAT  -  Nature.  Contour  interval  i  .s 
0.5  ms’  ,  negative  values  are  dashed. 
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Figure  7.22  (continued). 


LRTITIDE  (  BCCPEC3  ) 


LPTITUnE  C  OECMEES  ) 


Zonal  time  averaged  temperature.  (a)  Nature,  (b)  STATSAT , 

(c)  STATSAT  •  Nature,  (d)  WINDSAT,  (e)  WINDSAT  -  Nature, 

(f)  SSMSAT,  (g)  SSMSAT  -  Nature.  In  analyses,  values  (in  K)  li.ive 
been  multiplied  by  0.1,  contour  interval  0.5  (K/10) .  In 
differences  values  are  unsealed,  contour  interval  is  IK. 
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Fipure  7.24  Zonal  time  averaged  relative  humidity.  (a)  Nature,  (b)  STATSAl 
(c)  STATSAT  -  Nature,  (d)  WINDSAT,  (e)  WINDSAT  -  Nature, 

(f)  SSMSAT,  (g)  SSMSAT  -  Nature.  Contour  interval  is  5  percni 

negative  values  are  dashed. 
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Figure  7.25  Forecast  rms  error  growth,  500  mb  height.  (a)  Northern 

hemisphere  extratropics,  (b)  Tropics,  (c)  Southern  hemisphere 
extratropics . 
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Figure  7.26  Forecast  rms  error  growth,  200  mb  wind  vector.  (a)  NorLhciii 
hemisphere  extratropics,  (c)  Southern  hemisphere  extratropics 
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Figure  7.27  Forecast  rms  error  growth,  850  mb  relative  humidity. 

(a)  Northern  hemisphere  extratropics,  (b)  Tropics,  (c)  Soutlu  rii 
hemisphere  extratropics. 
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Figure  7.28 


Forecast  rms  error  growth.  (a)  850  mb  relative  humidity,  global 
statistics,  (b)  850  mb  cloud  cover,  global  statistics. 
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Figure  7.29  Forecast  rms  error  growth,  50  mb  height,  northern  hemisphore 
extratropics . 
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Figure  7.30  Calibrated  500  mb  height  rms  errors  for  southern  hemisphere 
extratropics.  (a)  Winter  OSEs,  (b)  Summer  OSE's. 
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I  Calibrated  200  mb  wind  vector  rms  errors  for  tropics.  (a 
OSEs,  (b)  Summer  OSEs. 
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Figure  7.32  Calibrated  850  mb  relative  humidity  sununer  OSSE  rms  errors  for 
northern  hemisphere  extratropics.  (a)  Relative  humidicv, 

(b)  Cloud  cover. 
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Figure  7.33  Forecast  500  mb  temperature  bias  for  northern  hemisphere 
extratropics . 
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Figure  7.34  Forecast  850  mb  relative  humidity  bias  for  northern  hemi.sphore 
extratropics . 
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8  San’.irarv  and  Conclusions 

Ihc  impact  of  the  183  GHz  temperature  and  RH  retrievals  on  the  AFGL  GDAS 
was  evaluated  by  comparing  the  results  of  OSSEs  using  the  SSM  data  instead  of 

(SSMSAT)  or  in  addition  to  TOVS  data  (SSM+TOVS)  with  baseline  OSSEs  using  the 

full  FGGE  data  set  (STATSAT)  and.  for  calibration  purposes,  no  satellite  data 
a:  all  (f.’OSAT)  .  Our  principal  findings  are  that  the  RH  forecast  errors  were 
slightly  affected  in  the  Northern  Hemisphere  extratropics,  but  a  marked 
i;.:f  rovement  could  be  seen  in  the  tropics  and  the  Southern  Hemisphere 
extratropics.  Interestingly,  though,  UINDSAT  moisture  forecasts  were  superior 
to  SS>!SAT  in  the  extratropics ,  even  though  only  RAOB  moisture  data  were  u.sed 

in  WINDSAT,  reflecting  the  dominant  role  of  the  mass  and  wind  fields  in 

forcing  the  moisture  field  in  the  extratropics. 

Gridpoint  to  gridpoint  comparisons  of  the  OSSE  anal;y.ses  with  the 
cor  re  .spending  nature  data  showed  the  SSMSAT  height  analyses  to  be  superior  to 
STATSAT  and  SSM+TOVS,  when  measured  in  terms  of  global  rms  errors.  Moisture 
analyses  were  clearly  improved  by  the  addition  of  SSM  data;  rms  errors  of 
SSMSAT  and  SSM+TOVS ,  which  were  almost  identical,  were  approximately  5%  less 
than  those  of  STATSAT  and  NOSAT  (22%  vs.  27%).  SSM+TOVS  and  STATSAT  had 
ccir.parable  height  errors.  Weather  features  were  identified  by  subjective 
synoptic  evaluations  which  were  analyzed  better  in  SSMSAT  than  SSM+TOVS.  The 
.s.-^mewhat  surprising  result  that  adding  TOVS  data  led  to  a  degradation  of  the 
height  analyses  was  caused  by  a  combination  oi  factors,  some  of  which  were 
related  to  the  anchoring  of  satellite  thicknesses,  others  to  data  selection 
and  quality  control  procedures. 

Cloud  cover  estimates  derived  from  the  relative  humidity  fields  were  too 
high,  Either  the  model  is  too  moist  or  the  relative  humidity  to  cloud  cover 
algorithm  needs  to  be  tuned.  In  any  case,  cloud  cover  differences  or 
comparisons  are  still  useful  since  all  relative  humidity  fields  converted  to 
cloud  cover  will  be  too  cloudy  in  the  same  way.  the  comparisons  of  rms 
differences  of  cloudiness  yield  the  same  results  as  comparisons  of  rms 
differences  of  relative  humidity.  This  might  have  been  expected  in  view  of 
the  facts  that  the  statistics  calculated  the  averages  over  large  samples  and 
that  the  relative  humidity  to  cloud  cover  relationship  is  simple.  This 
lelationship  is  nonlinear  but  it  is  1-to-l  and  monotonic  and  it  does  not 
depend  on  any  other  model  parameters  known  to  impact  cloudiness  such  as 
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temperature  lapse  rates,  vertical  wind  shears,  vertical  velocity  or 
divergence . 

The  OSSE  analyses  were  all  found  to  be  noisier  than  the  corresponding 
nature  data.  We  also  found  that  the  the  AFGL  model  has  a  tendency  to  warm  and 
dry  out  relative  to  the  ECMWF  nature  model.  This  warming  was  seen  in  all  the 
forecasts.  We  note  that  the  version  of  the  AFGL  model  which  we  used  has  no 
radiation  parameterization  and  hence  no  cooling  mechanism  although  there  is  a 
constant  source  of  warming  due  to  the  fixed  sea  surface  temperature. 

To  calibrate  these  OSSE  results  with  previously  conducted  OSEs,  error 
statistics  were  also  computed  between  the  OSSE  forecasts  and  the  simulated 
radiosonde  observations .  Impacts  in  the  N.H.  forecasts  were  larger  when 
calculated  as  grid  point  rms  errors  than  when  calculated  as  radiosonde  rms 
difference.  As  seen  in  the  synoptic  charts  the  greatest  impacts  tend  to  be 
polar,  however  the  verifying  radiosondes  tend  to  be  midlatitude.  The  impact 
of  the  SSM  data  was  found  to  be  generally  small  in  the  Northern  Hemisphere, 
but  quite  substantial  in  the  Southern  Hemisphere.  This  result  is  in  agreement 
with  numerous  previously  conducted  OSSEs  and  OSEs  which  measured  the  impact  of 
satellite  data.  When  measured  in  terms  of  500  mb  height.  Southern  Hemisphere 
forecast  errors  of  SSMSAT  (SSM+TOVS)  lagged  those  of  STATSAT  by  about  1'  (8) 
hours.  By  comparison,  addition  of  lidar  wind  profiler  data  to  STATSAT  led  to 
a  gain  of  36  hours  of  forecast  time.  These  forecast  time  gains  were 
essentially  unchanged  when  the  OSSE  results  were  calibrated  by  the  ratio  of 
OSSE  vs  OSE  STATSAT-NOSAT  impact.  The  RH  impact  could  not  be  calibrated  with 
OSE  results  because  there  is  essentially  no  difference  between  STATSAT  and 
NOSAT  OSE  moisture  errors. 

There  are  some  important  caveats  that  apply  to  the  results  reported 
here.  As  is  the  case  with  all  OSSEs  and  OSEs,  the  measured  impacts  apply  to 
the  particular  data  assimilation  system  used  here.  While  the  assimilation 
system  is  reasonably  "state  of  the  art",  some  aspects,  in  particular  the 
anchoring  of  the  satellite  thicknesses  and  the  limitations  of  data  selection, 
may  limit  the  extent  to  which  our  conclusions  are  generally  valid.  Since  we 
used  simulated  data  in  our  experiments,  the  realism  of  our  OSSE  results 
depends  on  how  realistic  the  simulated  observation  errors  were.  We  Cook 
particular  care  that  they  were  of  sufficient  size  and  had  the  appropriate 
error  correlation  structure.  Finally,  the  calibration  of  our  OSSE  results 
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wir!;  .'K  RAOB  statistics  suffprs  from  the  usual  problems  of  the  uneven 

.i is t r i hut i on  of  radiosondes,  particularly  the  bias  toward  land  areas  (where  no 

s  loll  ire  data  were  used),  and  the  small  sample  sizes  over  the  Southern 

hi-:r  ! phure  . 

Tirrre  is  cor;s i de table  opportunity  to  improve  and  refine  the  experiments 
reported  here  and  elsewhere.  Such  efforts  would  allow  the  quantification  od 
the  relative  impact  of  proposed  advanced  temperature  sounders  and  doppler  wind 
I.Lciars.  In  addition  cost  benefit  analyses  of  observational  accuracies  could 
ho  s’upported  by  such  studies.  It  is  also  i  rT’]i')r  t  ant  that  the  as.s  i  mi  1  at  ion 
sy.srcT  be  modified  to  best  take  advantage  of  (he  novel  observations.  And  in 
future  studies  it  will  be  important  to  carefully  simulate  the  geographical 
coverage  and  error  characteristics  of  proposed  instruments.  In  particular, 
(natural  pher.orr.ena  which  give  rise  to  correlated  observational  errors  should  be 
Included  to  the  extent  possible. 

Th.e  etror  chara>.:teriscics  and  distribution  of  simulated  temperature  and 
Rli  ’'ctrie.’als  for  example  should  depend  on  the  global  distribution  of 
clouds.  These  geophysical  parameters  are  in  turn  associated  (correlated)  in 
tin?  reel  world  with  the  geophysical  parameters  which  are  to  be  measured.  This 
is  oTie  cause  of  spatial  and  temporal  error  correlations  and  as  such  it.  should 
be  included  in  our  simulation  of  observational  errors.  Hence,  it  is  important 
th.at  the  model  used  to  generate  nature  also  can  provide  a  realistic 
description  of  the  cloud  field.  In  future  work  we  anticipate  using  a  newer 
nature  run  generated  by  the  ECMWF  using  T106  spectral  truncation  and  more 
coir:piete  physical  paramecerizations .  This  nature  run  includes  many  diagnostic 
fields  generated  by  the  physical  parameter izat ions . 

Ko  existing  global  model  has  fine  enough  resolution  to  represent  all 
scales  of  motion  which  exist  in  nature.  In  fact,  the  smallest  scales 
repiesented  by  models  are  usually  severely  damped  for  computational  reasons. 

In  real  data  assimilation  these  small  scales  are  considered  part  of  the 

oh;  e  rvat  lonal  error.  In  fact,  for  radiosondes,  this  is  the.  greatest  source  of 

error.  Although  the  length  scale  is  small,  these  errors  are  correlated. 

Th.c  T  (•  f  >0  re  simulated  observmt  ions  .should  include  spatially  corielatcd  errors 
Ti  r  ;  could  be  accomplished  by  unfiltering  the  nature  run,  thereby  restorinp. 
th-  (llest  resolvable  scales  as  .suggested  by  Hoffm.in  (  1  988). 
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Appendix  A.  MICROWAVE  SOUNDER  FIEIDS -OF- VI EU  FOR  THE  183  GHz  OSSE 


A.l  Approach 

A  basic  element  of  the  183  GHz  observing  system  simulation  experiment 
(OSSE)  is  the  inclusion  of  simulated  retrievals  of  water  vapor  and  temperature 
from  the  SSM/T-2  and  SSM/T-1  millimeter  moisture  and  microwave  temperature 
sounders,  respectively.  These  sensors  along  with  the  SSM/I  microwave  imager 
will  define  the  special  sensor  complement  relevant  to  meteorology  to  be 
carried  aboard  the  polar  orbiting  DMSP  satellites  in  the  1990's  time  frame. 

In  order  to  locate  the  additional  initialization  data  potentially  available 
from  this  observing  system  in  space  and  time  for  the  purpose  of  the  OSSE,  it 
is  necessary  to  develop  a  simple  sensor  scan  line  generation  program  (SLGP) 
for  these  sen;  rs .  The  approach  consists  of  a  simple  model  assuming  nominal 
characteristics  both  for  the  DMSP  orbital  platform  and  the  SSM/T-1  and  SSM/T-2 
cross -track  scanning  sensors.  Using  basic  information  regarding  both  platform 
and  sensor  and  latitude  and  longitude  of  each  possible  instrument  specific 
f ield-of -view  (FOV)  is  calculated  along  with  the  time  associated  with  the 
observation.  The  resulting  output,  therefore,  provides  a  simulated  data 
ephemeris  for  each  sensor  which  can  be  used  to  locate  the  available  satellite 
data  for  assimilation  into  the  numerical  weather  prediction  (NWP)  model.  Once 
the  potentially  available  data  is  located  in  this  manner,  the  appropriate 
errors  in  vertical  temperature  and  moisture  profiles  may  be  assigned  based  on 
their  dependence  on  such  factors  as  latitude  and  surface  type  (i.e.  ocean  vs. 
land)  . 

The  problem  is  simplified  considerably  by  assuming  that  the  DMSP  orbit  is 
circular,  the  earth  is  a  sphere,  and  the  satellite  subtrack  can  be  assumed  to 
be  a  great  circle.  In  this  case  concepts  from  spherical  trigonometry  can  be 
applied  to  both  subtrack  and  sensor  scan  line  location  calculations.  The 
calculation  should  thus  provide  a  useful  simulation  of  an  operational  sensor 
ephemeris.  The  calculation  may  be  divided  into  two  parts:  (a)  location  of 
the  satellite  subtrack  position  when  the  sensor  (i.e.  SSM/T-1  or  SSM/T-2)  is 
observing  a  specified  FOV,  and  (b)  location  of  the  center  of  the  FOV  itself. 
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A. 2  Platform  and  Sensor  Characteristics 


It  is  assumed  that  the  DMSP  platform  is  in  a  circular  orbit  at  an 
altitude  of  833  km  with  a  nominal  period  of  101.35  minutes.  The  Inclination 
of  the  orbital  plane  with  the  equatorial  plane  of  the  earth  is  98.7  degrees 
insuring  sun  synchronic ity.  Other  key  satellite  parameters  are  summarized  in 
Table  A-l.  These  include  the  time  and  longitude  of. the  ascending  node  chosen 
to  start  the  calculation.  The  ascending  node  positions  the  satellite  at  the 
equator  crossing  from  southern  to  northern  hemisphere.  The  longitude  and  time 
of  the  ascending  node  determine  the  local  solar  time  at  which  the  satellite 
transmits  the  meridian.  This  determines  whether  data  is  being  collected  from 
a  morning  or  noon  satellite. 

Both  the  SSM/T-1  and  SSM/T-2  sounders  are  cross -track  scanning  sensors. 
Data  coverage  is  provided  by  scanning  perpendicular  to  the  satellite  subtrack 
in  a  specifiable  scanning  pattern.  The  desired  scanning  pattern  consists  of  a 
fixed  number  of  scan  positions  incremented  along  the  scan  and  an  InterScan 
calibration  period.  Each  scan  position  determines  a  f ield-of-view.  It  should 
be  noted  that  modeled  FOV's  are  approximate.  One  source  of  error  are  those 
associated  with  pointing  angle  accuracy.  For  example,  the  beam  position 
repeatability,  scan  plane  alignment,  and  distribution  of  the  actual  he§m 
position  about  the  desired  beam  position  have  tolerances  of  about 
0.4  degrees.  Additionally,  the  beamwidth  for  the  i>SM/T-2  used  in  our 
calculation  corresponds  to  the  183  GHz  channels  only.  The  two  lower  frequency 
channels  at  91.6  and  150  GHz  have  beamwidths  of  3.7  and  6.7  degrees, 
respectively.  Finally,  since  beamwidths  for  both  sensors  are  determined  by 
the  3  dB  power  points,  there  is  some  overlap. 

The  SSM/T-1  views  7  positions,  one  at  nadir  and  three  to  either  side. 

The  scan  angle  increment  is  12  degrees  and  the  FOV  is  twelve  degrees  for  a 
total  scan  from  -36  to  36  degrees  of  nadir.  The  total  swath  width  is  thus 
42  degrees  to  either  side  of  nadir.  The  SSM/T-2  has  a  total  of  28  viewing 
positions  which  scan  from  -40.5  to  40.5  degrees  of  nadir  in  steps  of  3 
degrees.  Each  FOV  is  approximately  3  degrees  so  that  the  total  swath  width  is 
identical  to  that  of  the  SSM/T-1.  The  SSM/T-2  never  actually  views  the  nadir 
position  along  the  satellite  subtrack.  Two  scan  positions,  those  denoted  14 
and  15  view  -1.5  and  1.5  degrees,  respectively,  from  nadir.  The  scan  patterns 
for  the  SSM/T-1  and  SSM/T-2  sensors  are  illustrated  in  Figure  A-l. 
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Thf  total  scan  cycle  period,  T' ,  for  Che  SSM/T-1  sensor,  l.e.  the  time  it 
ta'fCi  s  to  sweep  on  a  full  360  degrees  including  FOV  sampling  and  calibration 
is  32  seconds.  Tlie  scan  cycle  period  for  tlie  SSM/T-2  is  8  seconds.  During 
rhf  scan  cycle  period  each  sensor  views  FOV's  for  42  degrees  and  calibrates 
etc.  for  the  rest  of  the  cycle.  Ass'uming  a  constant  rotation  rate  for  the 
Cities,  a  FOV  sampling  time  and  an  interscan  period  can  ho  calculated.  These 
and  other  sensor  related  characteristics  are  given  in  Table  A-2.  In  ortier  to 
keep  account  of  the  sensor  position  at  any  time,  we  define  a  scan  line  in.dex, 
mJ  ,  and  a  scan  position  index,  n^ ,  for  each  sensor  j ,  where  j  =  1,2  denote  the 
SSM/T-l  and  SSM/T-2,  respectively.  Associated  with  each  index  pair  (rn^  ,  n-^  ) 


will  h« 


'.gitude,  latitude,  and  time  (A,(^,t).  The  calculation  is 


initialized  by  specifying  the  time  and  longitude  oi.  the  ascending  node  and  the 
ir.icial  sensor  positions.  (By  definition,  the  latitude  at  the  ascending  node 
tiite  is  Zero.)  Tlie  initial  sensor  positions  are: 


1  1 
V  i.t  .  n 


( m  ,  n  )  C  ( 


1  1 
('"0  ,  IIq  ) 


2  2 


(A.  1 ) 


ar.d  provide  a  starting  point  for  simulating  the  scan  cycle.  The  scan  ^)Osition 
index  is  incremented  by  one  until  the  total  number  of  scan  positions,  N-^  ,  for 
each  sensor  is  reached.  At  this  point,  the  scan  line  index  is  incremented. 
Based  on  these  indices  the  sensor  scan  angle  and  time  past  the  ascending  node 
can  be  evaluated.  The  scan  angle  at  scan  position  is  given  by: 


+ 
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nJ-INT(^  )  +  1 
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(A. 2) 


where  (/J  is  the  initial  scan  angle,  n-^  is  the  scan  position  index,  is  the 
total  number  of  scan  positions,  and  is  the  scan  angle  increment  for 

sensor  j.  Scan  angles  can  take  on  values  which  are  both  positive  and  negatiw 
depending  whether  they  are  to  the  starboard  or  port  side,  respectively,  of  rhe 
spacecraft  velocity  vector.  The  time  along  the  satellite  .subtrack  can  be 
evaluated  by  noting  the  elapsed  time  since  the  ascending  node  indexed  by  the 
number  of  scan  lines  acquired  and  the  current  scan  position.  Thus  the  current 
time  will  be : 


(A. 3) 


t(m^,n-^)  -  tg  +  { (N^  -  n^)  +  n^](t^/N-^)  +  -  m|)  At^ 

+  (m^  -  mg  -  DN"^ 


where  is  the  FOV  sampling  time  and  A-^  is  the  InterScan  period. 


Table  A-1.  DMSP  Satellite  Platform  Characteristics 


Parameter 

Description 

Assumed  Value 

P 

Satellite  period  (m) 

101.35  min 

€ 

Satellite  inclination  to  equatorial 
plane  (deg) 

98.7  deg 

Zq 

Satellite  altitude  (km) 

833  km 

Rq 

Radius  of  the  Earth  (km) 

6371  km. 

^0 

Time  of  ascending  node  (h.m.s) 

0:0:0  (GMT) 

^0 

Longitude  of  the  ascending  node  (deg) 

0.0  deg 

1  Jib  1  t\  ■  . 


St-nsor  Characteristics  for  the  SSM/T-1  and  SSM/T-2  Sensors 


P.i  r-ii^e  t  e  r 

Description 

Assumed 
SSM/T- 1 

Val  uc 
SS.M/T  -  : 

O 

Initial  scan  angle  sensor  j  (deg) 

O 

o 

1  .  5 

b.u-j 

Cross -track  scan  angle  increment:  for 
sensor  j  (deg) 

12.0 

3.0 

J 

Total  number  of  scan  positions  for  sensor  j 

7 

28 

T 

Scan  period  for  sensor  j  (s) 

32.0 

8.0 

t-i 

FOV  sampling  time  for  sensor  j  (s) 

7.47 

1.87 

j 

InterScan  period  for  sensor  j  (,s) 

24.  b 3 

6.  13 

IT.-' 

Scan  1 ine  index  for  sensor  j 

Scan  position  index  for  sensor  j 

( nJ  ) 

scan  angle  at  scan  position  index  n^ 
for  sensor  j  (deg) 

ej  (riJ  ) 

Azimuth  angle  at  scan  position  index  nJ 
for  sensor  j  (deg) 

- 

-5 


Orbital  Track 


_ SSM/T-1 

- SSM/T-2 


Figure  A-1.  SSM/T-1  and  SSM/T-2  FOV  scan  line  pattern. 
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A.  :  r i o 1 d ' o f ■ V i cw  Location 

The  r.f'ometry  rentiiied  to  calculate  locations  for  each  f ield-of-view  is 
illustrated  iit  Figure  h-2.  Here  the  ascending  node  is  denoted  by  tha  letter 
A.  a  general  point  -ilong  the  satellite  subtrack  as  letter  B,  and  the  starboard 
tx^ent  of  a  ger.oral  scan  line  denoted  by  the  letter  E.  The  pole  is 
lctt-‘i  C.  Assuming  t  tie  satellite  subtrack  (arc  AB)  to  be  a  great  circle,  we 
r.:!';  c.ilrulate  tlie  latitude  and  longitude  of  point  B  bv  applying  Napier's 
,in,. logics  i,Selhv,  l^itVt,  p.  2,^3)  for  the  spherical  triang,lc  ARC.  The  latitudi- 
at  point  B  will  be  given  by  the  arc  length  (in  dc gt •  ■  s)  of  side  a: 

\(t)  -  QO.n  -  (A.  4) 

wi.il  e  the  longitude  will  be  given  bv  the  angle  C: 

0 '  :  -  d,,  t  C  +  ("  t  -  c w  (A.  ii) 

o  U 

where  u  is  the  angular  velocity  of  the  earth.  The  longitude  increment  is 
addetl  to  that  of  the  ascending  node  and  a  correct:ion  is  made  for  the  rotation 
of  the  earth  since  the  ascending  node  time. 

The  unknowns  a  and  C  are  obtained  from  our  knowledge  of  sides  c  and  b  and 
angle  .A.  Angle  A  is  a  con.stant  and  is  equal  to  the  differenct-  in  angle 
between  the  orbital  plane  of  the  satellite  and  a  typical  meridian  or: 

A  =■  e  -  90.0  -  8.7  deg  (A. 6) 

Side  b  is  simply  the  arc  length  from  point  A  (on  the  equator)  to  the  pole 
point  C,  or: 


b  -  90.0  deg  (A. 7) 

The  length  of  side  c  is  obtained  from  the  elapsed  time  since  the  ascending 
node  realizing  that  one  orbital  period  corresponds  to  360  degrees.  Thus: 

c  -  360  (t  -  tQ)  /  p  (A. 8) 

a  -  2  tan'^  [h(b,c,B,C)l  (A. 9) 
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C  “  Can'^  [f(b.c,A)]  -  tan'^  [g(b,c.A)] 


(A. 10) 


where 


B  -  tan'^  [f(b,c,A)]  -  tan'^  [g(b,c,A)] 
and  the  functions  f,  g  and  h  are  given  by: 


(A, 11) 


f(b,c,A)  -  tan  ^  (B+C)  -  cos  ^  (b-c)  sec  ^  (b+c)  cot  A  (A. 12) 

g(b,c,A)  -  tan  ^  (B-C)  -  sin  |  (b-c)  esc  ^  (b-i-c)  cot  ^  A  A. 13) 

h(b,c,B,C)  -=  tan  ^  a  “  Can  ^  (b-c)  sin  ^  (B-hC)  esc  ^  (B-C)  (A. 14) 

Upon  substitution  of  the  quantities  A  (A. 6),  b  (A. 7),  and  c  (A. 8)  into 
equations  (A. 11)  through  (A. 14),  the  unknowns  a  (A. 9)  and  C  (A. 10)  can  be 
evaluated.  Equations  (A. 4)  and  (A. 5)  then  yield  the  desired  subtrack  point 
latitude  and  longitude,  respectively. 

Similar  considerations  apply  to  the  location  of  specified  f ields - of - view 
along  the  scan  line.  In  this  case  the  spherical  triangle  BEC  (Figure  A-2)  is 
considered.  The  length  of  .side  a  (equation  A. 4)  is  simply: 

a  -  90.0  -  A(t)  (A. 15) 

The  angle  at  B  is  determined  from  the  orbital  inclination  and  is; 

B  -  81 . 3  deg  (A. 16) 

The  length  of  the  scan  line  segment  f^J  at  position  n^  in  degrees  can  be 
obtained  from  the  scan  angle  ('A.2)  and  the  local  zenith  angle  given  by; 


A-8 


_  i  .  -I  %  ^  h 

0  -  Sin  - ;; - 

n  Rq 


sin  "t 


(A.W) 


From  Figure  A- 3,  it  can  be  seen  that  this  angle  is: 


f  J  «  e  J  -  >1/ 
n  n  n 


(A. 18; 


The  Napier  analogies  (equations  (A. 9)  through  (A. 14))  can  now  be  applied  to 
solve  for  the-  latitude  and  longitude  for  the  FOV  specified  by  each  sensor's 
scan  indices.  The  unknowns  are  the  length  (in  degrees)  of  the  side  d  and  tin 
new  angle  C  included  by  ECD  in  Figure  A-2.  Given  these,  the  latitude  and 
lonpitude  of  tTie  f ield-of-view  are: 


A  ( n"'  ,  ) 


(A. 19) 


i  ( n'^  .  ) 


(A. 20) 
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Figure  A-4.  Subsatellite  trade  for  polar  orbiting  platform  (24  hours  of 
orbits) . 
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